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INTRODUCTION 



Theories for the form ation and e volution of planets depended primarily on geophysical data 
from the Solar System (see BrushlEoQQ . and references therein). Starting in the 1940's, astrophysical 



data began to provide new insights. Discoveries of pre-main sequence stars in Taurus-Auriga, 
Orion, and other reg ions led to the concept that stars form in giant clouds of gas and dust (see 



r^: 

Kenvon et al.ll2008bl . and references therein). Because nearly every young star has a circumstellar 



disk with enough mass to make a planetary system, theorists began to connect the birth of stars 
to the birth of planets. Still, the Solar System remained unique until the 1990's, when the first 
discoveries of exoplanets began to test the notion that planetary systems are common. With 
thousands of (candidate) planetary systems known today, we are starting to have enough examples 
to develop a complete theory for the origin of the Earth and other planets. 

Here, we consider the physical processes that transform a protostellar disk into a planetary 
syste m around a single star. Instea d of discussing the astrophysical and geophysical data in detail 



(e.g.. iDauphas &: Chaussidonll201ll ). we focus on a basic introduction to the physical steps involved 
in building a planet. To begin, we discuss several observational constraints from the wealth of 
astrophysical and geophysical material in ^ We then describe the global physical properties and 
evolution of the disk in ^ Aside from the special conditions required for fragments in the disk 
to collapse directly into giant planets ( §5.3p . most planets probably grow from micron-sized dust 
grains. Thus, we consider how a turbulent sea of grains produces the building blocks of planets, 
planetesimals (® , and how ensembles of planetesimals collide and merge into planets which may 
later accrete a gaseous atmosphere (^. We conclude with a brief summary in ^ 



2. OBSERVATIONAL CONSTRAINTS ON PLANET FORMATION THEORIES 

We begin with the main observational constraints on planet formation processes, including raw 
materials, timescales, and outcomes. Detailed studies of the Solar System and the disks around the 
youngest stars yield strict limits on the mass available and the time required to make a planetary 
system. The diverse population of exoplanets illustrate the many outcomes of planet formation. 



2.1. Lessons from the Solar System 

Until the discovery of exoplanets, the Solar System was the only known planetary system. The 
Solar System will continue to provide the most detailed data on planet formation, despite obvious 
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issues of statistical significance and anthropic bias. 



2.1.1. The Solar Nebula 



The alignment of major planets in the ecliptic plane suggests that they formed within a flat- 
tened disk or "nebula." The philosopher Immanuel Kant and the mathematician Pierre-Simon 
Laplace are often credited for this "nebular hypothesis." However, the scientist and theologian 
Emanuel Swedenbourg first recorded this insight in his 1734 Principia. For a long time the nebular 
hypothesis competed with the theory, proposed by the naturalist Buffon, that planets were tidally 
extracted from the Sun during an encounter. Though Laplace dismissed the encounter theory as 
being inconsistent with the circular orbits of the planets, it su rvived to reac h peak popularity in 
the early 20th century as the Chamberlin-Moulton hypothesis ( Jeffreys 19291 ). Once Payne ( 19251 . 
1925b) identified hydrogen as the most abundant element in stars, the nebular hypothesis regained 
favor. Adding the tidal theory's idea of planetesimals - small solid particles that condense out of 
hot gas ~ the nebular hypothesis began to develop into a robust theory for planet formation. 

The minimum mass solar nebula (MMSN) provides a simple estimate of the mass available in 
planet forming disks. The recipe for the MMSN is to distribute the mass currently in the solar 
system's planets int o abutting ann uli, adding vo l atile elements (mainly hydrogen gas) until the 
composition is Solar. Kuiper ( 195ll ) and Cameron ( 1962 . se e his T able 5) estima t ed the mass of the 
MMSN as a few percent of a solar mass. Weidenschilling ( 1977bl ) and Hayashi ( 1981 ) fit the now 
canonical surface density law, bravely smoothing the mass deficits in the regions of Mercury, 

Mars, and the asteroid belt, and the abundance uncertainties for the giant planets. The roughness 
of the fit is immaterial: the MMSN is not a precise initial condition, but a convenient fiducial for 
comparing disk models. 

We use the same MMSN as IChiang YoudinI (|20ld l with disk surface density profiles: 
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where subscripts "g" and "p" respectively denote gas and particles (condensed solids) and R is 
the distance from the Sun. The parameter F scales the total mass; F = 1 is a reference MMSN. 
Integrated out to 100 AU, the MMSN disk mass is O.O3M0. 

The parameter Z^ei scales the ratio of solids-to-gas, or disk metallicity as 



■^disk = ^p/^g = O.OlSZrel . 



(2-3) 



evolves during the planet formation process, as evidenced by the enrich ment of heavy elements 
relative to H in Jupiter and Saturn. Our fiducial value is normalized to the Lodders ( 20031 ) analysis 
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of (proto-)Solar abundances, which can be approximately summarized as: 



1 T < 40 K 

0.78 40K<r<180K 

0.33 180 K < T < 1300 — 2000 K 

r> 1300 — 2000 K 



(2-4) 



The abundance of soHds decreases with increasing temperature due to the subhmation of (most 
significantly) methane ice above 40 K, water ice above 180 K, and dust over a range of temper- 
atures from roughly 1300 — 2000 K, covering the condensation temperatures for different min- 
erals. Note that this definition of disk metallicity ignores the heavy elements in the gas phase, 
which are at least temporarily unavailable to produce planetary cores. For the disk temperature, 
this work adopts eg. (|3-29p . the result for an irradiated disk with a self-consistently fiared surface 



(jChiang &: Goldreich 



19971 1 



In the Solar System today, the temperatures of solids are set by an equilibrium between heat- 
ing and cooling. Although gravitational contraction (Jupiter) and tidal heating (many satellites, 
including the Moon) contribute some heating in a few objects, the Sun is the primary source of 
heating. Radiation from the Sun peaks at a wavelength of Xmax,Q ~ 0.5 /um. Objects with radius 
r, peak wavelength Xmax,g < f-, and no atmosphere radiate as nearly perfect blackbodies. Equating 
the energy they receive from a central star (7rr^L*/47ri?^) with the energy they emit (47rr^(7sBT'gg j,^) 
leads to an equilibrium temperature , 
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where Lq is the luminosity of the Sun. Small grains with r > 1 /im and r < Xmax,g 
inefficiently. In most cases, the radiative efficiency is 
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where Aq is the critical wavelength and q ^ 1-2 depends on grain properties. Usually Aq is roughly 
equal to the grain radius r. Because they can only radiate efficiently at short wavelengths, these 
grains have much larger temperatures. For q = 1, 
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To derive eq. (j2-7p . include the efficiency in the grains' emitted radiation (oc 1/e) and relate the 
wavelength of peak emission to the grain temperature with Wein's Law (A oc l/Tgg^s). 

Coupled with the condensation temperatures in eq. (|2-4p . these definitions allows us to identify 
the "snow line." Also known as the "frost line," this annulus in the dislfl which separat es an 
inner region of rocky objects from an outer region of icy objects ( Kennedv &: Kenvon 20081 ). For 



or, more generally, a spherical shell surrounding the central star 
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blackbody grctins, tli6 Welter condcnscition temperature of 180 K implies Rsnow 

K, 2.7 AU, roughly 

coincident with the asteroid belt. Similarly, the methane condensation temperature of 40 K yields 
another region beyond the outer edge of the Kuiper belt at 48 AU where solid objects have a 
combination of water and methane ice. Inside of ~ 0.1 AU, dust evaporates; rocky grains cannot 
exist so close to the Sun. 



2.1.2. Isotopic Timescales 



Meteorites delivered to Earth from the asteroid belt provide the most detailed chronology of 
the early Solar System (e.g., bauohas fc Chaussido-^ B. Primitive meteorites from asteroids 
that did not undergo differentiation (or other significant alteration) preserve the best record of their 
formation. These primitive meteorites are called chondrites because they contain many chondrules. 
Chondrules are glassy inclusions, with a typical size ~ 0.1 — 1mm. They provide evidence for high 
temperature melti ng events in the Solar nebula. The nature of these melting events is debated and 
beyond our scope ( Connolly et al. 20061 ). Calcium Aluminum Inclusions (CAIs) are also present in 
primitive meteorites. CAIs experienced even more extreme heating than chondrules. 

200^), CAIs are the oldest known objects 



With ages up to 4567.11 ± 0.16 Myr (|R.ussen et al 



in th e Solar System. This age is consistent with current results for the main sequence age of the 



Sun (Bonanno et al 



2003). The absolute ages of CAIs are measured by lead-lead dating, which 
makes use of half-lives, ii/2) of uranium isotopes that are conveniently long. The decay chain of 
235u ^207p]3 i^g^ ^^^^ ^ 0.704 Gyr, while ^ssu ^206p|3 j^^^g ^^^^ ^ 4 47 q^j. 

Radioisotopes with short half-lives yield accurate relative ages of meteorite components. These 
isotopes are "extinct;" they have decayed completely to daughter products whose abundances rela- 
tive to other isotopes result in an age. The extinct isotope ^^Al decays to ^^Mg in ti/2 = 0.73 Myr. 
The abundance of ^^Mg relative to ^^Al and to ^^Mg yields an age for the sample. Aside from its 
use as a chronometer, ^^Al is a powerful heat source in young protoplanets. 

Both absolute (lead-lead) and relative (^^Al) ages support a planet formation timescale of a 
few Myr. Most CAIs fo rmed in a narrow w indow of 1 — 3 x 10^ yr; chondrule formation persisted 
for ~ 4 Myr or longer. iRnssell et all tood ) discuss systematic uncertainties. Here, we emphasize 
the remarkable agreement of the few Myr formation time derived from primitive meteorite analyses 
and protoplanetary disk observations ( §2.2p . 

The assembly of terrestrial planets from planetesimals requires tens of Myr. Isotopic analysis 
of differentiated Solar System bodies (including Earth, Mars and meteorites) probe this longer 
timescale. The decay of radioactive hafnium into tungsten, ^^^Hf — j-^^^ with ti/2 = 9-8 Myr, 
dates core-mantle segregation. Tungsten is a siderophile (prefers associating with metals) while 
hafnium is a lithophile (prefers the rocky mantle); thus, Hf-W isotope ratios are the primary tool 
to date differentiation (Chapter by Barlow). Studies of Hf-W systematics indicate that asteroid 
accretion con tinued for ^ 10 My r, Mars' core formed within 20 Myr, and the Earth's core grew over 
30-100 Myr ( Kleine et al.l 20091 ). Astronomical observations of debris disks ( §2.2p and dynamical 



Astronomical observations of debris disks ( §2.2 
studies of terrestrial planet accretion (^ support these longer timescales. 
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2.1.3. Water 

After hydrogen, water is the most abundant molecule in disks (e.g., Najita et al. 200?! ). Inside 
the snow line, water exists in the gaseous phase, though it dissociates at T > 2500 K. As water 
vapor interior to the ice line diffuses past the snow line, it condenses into icy grains. The snow line 
thus acts as a cold trap, where the enhanced mass in water ice fea. | 2-4p. should accelerate the 
growth of planetesimals and perhaps gas giant planets (jStevenson h Lunindll988. ). 

Water is abundant throughout the Solar System (e.g., iRivkin et al.l 120021 : iJewitt et al.l 120071 . 
and references therein). Outside of the Earth, water appears in spectra of comets, Kuiper belt 
objects, and satellites of giant planets (including the Moon) and bound within minerals on Mars, 
Europa, and some asteroids. Because they can be analyzed in great detail, meteorites from asteroids 
provide a wealth of information on water in the inner Solar System. Many meteorites show evidence 
for aqueous alteration prior to falling onto the Earth. Most groups of carbonaceous chondrites 
and some type-3 ordinary chondrites contain hydrated minerals, suggesting association with liquid 
water. Despite some evidence that grains might react with water prior to their incorporation into 
larger solids, most analyses of t he mineralogy suggest hydr ation on scales of mm to cm within 
chondrites and other meteorites ( Zolenskv &: McSween 19881 ). 



The water content of Solar System bodies helps trace the evolution of the snow line (iKennedy &: Kenyon 



20081 ). Hydration within meteorite samples dem onstrates that the snow line was at least as close as 



2.5-3 AU during the formation of the asteroids (jRivkin et al.ll2002. ). Radiometric analyses suggest 
hydration dates from 5-10 Myr after the formation of the Sun, close to and perhaps slightly after 
the formation of chondrules. Ice on the Earth and on Mars suggests the possibility that the snow 
line might have been as close as 1 AU to the proto-Sun, a real possibility for passively irradiated 
disks ( ^3.4p . Within the terrestrial zone, the rise in water abundance from Venus (fairly dry) to 
Earth (wetter) to Mars and the asteroids (wetter still) points to processes that either distributed 
water throughout the inner Solar System (with a preference for regions near the snow line) or 
inhibited accretion of water (either vapor or ice) from the local environment. Abundance analyses, 
including D/H and noble gases, help to probe the (still imperfectly known) history of water in the 
inner Solar System. 



2.2. Disks Surrounding the Youngest Stars 

Observations of young stars provide additional constraints on the early evolution of planetary 
systems. In nearby star-forming regions (Orion, Ta urus, etc), nearly every star with an age of 1 
Myr or less has an optically thick circumstellar disk ( Williams &: Cieza 201 ll ). The disk frequency 



seems independent of stellar mass. The data suggest the disks are geometrically thin, with a vertical 
extent of roughly 10% to 20% of their outer radius. They are composed of molecular gas and dust 
grains with sizes ranging from a few microns up to several mm. Around solar-type stars, young 
disks have typical luminosities of order Lq and radii of order 100 AU. 

Estimating disk masses is challenging. Aside from a few transitions possible only in warm 
material near the central star, H2 is undetectable. The next most abundant molecule, CO, is 
optically thick and provides a crude lower limit to the total mass. Current estimates rely on 
a conversion from the dust emission at mm wavelengths to a dust mass and then to a mass in 
gas. These estimates are highly uncertain due to ignorance of dust-to-gas ratios and grain size 
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distributions. The typical assumptions give disk masses ~ 0.01 Mq, a factor of 2-4 smaller than 
the MMSN. The dispersion f or ensembles of a few hu ndred systems is roughly an order of magnitude 
( Andrews Sz Williams! 2005 : Williams Cieza l201lh . The size of the typical disk is similar to the 
semimajor axes of orbits in the Kuiper belt. 

Current data demonstrate that more massive young stars have more massive disks. Recent 
observational programs concentrate on whether the ratio of disk mass to stellar mass is roughly 
constant or increases with stell ar mass. Despite the larg er scatter in this ratio at every stellar mass, 



the relation is probably linear ([Williams &: Ciezal 12011 



High resolution radio observations reveal interesting limits on the distributions of surface den- 
sity and temperature within the brightest and most massive disks around nearby stars. Although 
disks with a broad range of surface density gradients ar e observed (U oc with n ~ —0.6 

to 1.5), most observations iii dicate a typical n ^ 0.5-1.5 ([Andrews &: Williamsl l2005l : llsella et al 



20091 : IWilhams fc Ciezall201lh . Thus the surface density gradient in the MMSN is steeper than the 



average protostellar disk but within the range observed in disks around other young stars. 

The evolution of protostellar disks sets severe limits on the timescale for planet formation. 
Optical and ultraviolet spectra of young stars show that material from the disk flows onto the 
central star. The rate of this flow, the mass accretion rate, drops from well in exc ess of 10~ 



Mq yr ^ at 1 Myr to much less than 10 Mq yr ^ at 10 Myr (jWilliams Sz Ciezall201lh . Declining 



mass accretion rates imply much less gas near the young star. At the same time, the fraction of 
young stars with opaque dusty disks declines from nearly 100% to less than 1%. Fewer dusty disks 
implies the solid material has been incorporated into large (km-sized or larger) objects, accreted 
by the central star, or driven out of the system by radiation pressure or a stellar wind. Direct 
constraints on the amount of gas left in 10 Myr old systems without opaque disks are limited to a 
few systems. 

Among older stars with ages of 3-10 Myr, many disks have substantial mass beyond 30-50 AU 
but have inner holes apparently devoid of much gas or dust. The frequency of these "transition" 
disks suggests the evolution opaq ue disk — )• opaque disk with mner hole — )• no opaq ue disk takes 
from 0.1-0.3 Myr up to 1-2 Myr (jCurrie fc Siciha-Aguilailboill : lEspaillat et ahlboid ). 

Once the opaque di sk disappear s, many pre- main and main sequence stars remain surrounded 
by 1-10 //m dust grains ( Wyatt 20081 ). This material lies in a belt with radial extent 6R ~ 0.1— 0.5i? 
and vertical height 6z ~ 0.05 — O.li?, with R ~ 1-100 AU. Infrared spectroscopy suggests the grains 
have compositions similar to material in the comets and the dust (Zodiacal Light) of the Solar 
System. The total mass, a few lunar masses, exceeds the mass of dust in the inner solar system 
by factors of 100-1000. These properties are independent of stellar metallicity and many other 
properties of the central star. However, the frequency of and the amount of dusty material in these 
dusty disks peaks for sta rs with ages of 10-20 Myr and then declines approximately inversely with 
time (jCurrie et al.ll2008l ). 

These disks place interesting limits on the reservoir of large objects around stars with ages of 
10-20 Myr. Among several possible grain remo val mechanisms, the most likely are radiation pro- 
cesses and collisions (iBackman &: Parescdll993l ). Grains orbiting the star feel a headw ind from the 



incom ing radiation from the central star, which causes the grain to spiral into the star (iBurns et al 



19791 ). If the grains have a mass density p,, the orbital decay time for this Poynting- Robertson 
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drag is 



t 



f 47rrp,\ / c^R"^ 



pr 



710 p. 



yr 



(2-8) 



3 / V 

The decay time, tpr ~ 1 Myr, is much shorter than the 100 Myr to 10 Gyr main sequence hfe- 
time of the central star. The large frequency of dusty disks among 50-500 Myr (1-10 Gyr) old 
A-type (G-type) stars suggests the grains are continually replenished over the main sequence life- 
time. By analogy wi th the Solar System, where trails of dust result from collisions of asteroids 
( Nesvornv et al.l 20031 ). high velocity collisions among large (but undetectable) objects can replen- 
ish the dust. Adopting typical sizes for asteroids, ~ 1-10 km, a reservoir containing ~ 10 Mq of 
material can explain the amount of dust around young stars and the time evolution of the dust 
emission among older stars. Because this mass is between the initial mass of solids in protostellar 
disks (~ 100-1000 M®) an d the dust in ass in the Solar System (< 10""^ M®), these systems are 
often called "debris disks" ( WyattI liooi . also Chapter by Moro-Martin) . 



2.3. The Exoplanet Revolution 



The current pace of exoplanet discovery is extraordinarily rapid. Here, we highlight the main 
insights exoplanets bring to theories of planet formation. 

Within the d iscovery space accessible with current techniques, exoplanets fill nearly all available 



phase space (e.g.. lCumming et al.ll2008l : iGould et al.ll2010l : iHoward et al.ll201ll : I Johnson et al.ll2011 

and references therein). Among ~ 10% to 30% of middle-age solar- type stars, exoplanets lie within 
a few stellar radii from the central star out to several AU. Because ther e are rnany ni ultiplanet 
systems, the sample of short period planets implies more planets than stars (|Youdinll2011bl ^. Though 
detection becomes more difficult, the frequency of exoplanets seems to grow with increasing distance 
from the parent star. The orbits have a broad range of eccentricities, with a peak at e ~ 0.2. Planet 
masses range from a rough upper limit at 10-20 Mj to a few Earth masses. Around stars with the 
same mass, lower mass planets are more frequent than higher mass planets. More massive stars 
tend to have more massive planets. 



There is some e vidence that exoplanets are more likely around metal-rich stars (lGonzaleall997l : 



Johnson et al 



20inl ). With the large samples available, this "planet-metallicity correlation" is now 
unambiguous for gas giants with masses ranging from the mass of Saturn up to ~ 10 Mj. Among 
lower mass planets, small samples currently prevent identifying a clear correlation. Larger samples 
with the Kepler satellite will yield a better test of the planet-metallicity correlation as a function 
of planet mass. 

The origin of any planet-metallicity correlation establishes some constraints on formation the- 
ories. If the initial metallicity of the disk is identical to the current metallicity of the star, then 
gas giants - and perhaps other planets - form more frequently in more metal-rich disks. However, 
planets could pollute the stellar atmosphere after the rest of the gaseous disk disperses, raising the 
metallicity of the star above the initial metallicity of the disk. In this case enhanced metallicity 
would be a result riot a cause of planet formati on. Current data contradicts the pollution hypothesis 
(jFischer &: Valentill2005l : iPasquini et al.ll2007l ). but more study of the diverse exoplanet population 
is warranted. 

To quantify the planet-metallicity correlation, Johnson et al. ( 20ld ) fit the frequency, /, of 
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giant planets as a joint power-law in stellar mass, M^, and metallicity, Zpe oc log;^o(Pe/H]), 

/ oc M^zi . (2-9) 

For giant planets in the California Planet Survey, a = 1.0 it 0.2 and (3 = 1.2 it 0.2. Ignoring 
the dependence with stellar mass (a = 0) introduces a bias, but yields a stronger relation with 
metallicity, (/3 = 1.7 it 0.3). For more massive stars with M^, > IAMq, a = 1.5 it 0.4 and 
P = 0.73 lb 0.35. Thus, the formation of giant planets around more massive stars is more sensitive 
to stellar mass and less sensitive to metallicity than for lower mass stars. 

With exoplanet samples growing so rapidly, new analyses will change at least some of these 
conclusions. The most firm conclusions - that exoplanetary systems are common and have nearly 
as much diversity as possible - provides a good counterpoint to the wealth of detail available from 
in situ analyses of the Solar System. 



3. DISK PROPERTIES AND EVOLUTION 



Stars form within collapsing clouds of gas and dust. When a cloud collapses, most infalling 
material has too much angular momentum t o fall directly onto the nascent protostar. This ga s forms 
a rotationally supported circumstellar disk ( Cassen &: Moosmanlll98ll : iTerebev et al.lll984l ). If the 
angular momentum in the disk is transported ra dially outward, gas can accrete on to the central 
star. Although jets launched near the protostar (|Shu et al.ll2"o0ol ) or from the disk (IPudritz et al 



20071 ) or both can remove significant angular momentum from the disk, most analyses concentrate 
on how angular momentum flows through the disk. 

At early times, the disk mass Mjisk is similar to the stellar mass M^. For a circumstellar disk 
with surface density X! and radial flow velocity vr, the rate material flows through the disk as a 
function of radial coordinate R is 

M = -IttRvrU . (3-1) 

Positive M corresponds to gas flowing towards, and eventually draining onto, the central star. If 
the mass infall rate from the surrounding molecular cloud Mi exceeds M, the disk mass grows. If 
Afnisk exceeds ^ O.SM^, gravitational inst abilities within the disk can produce a binary companion 
(| Adams et al.lll989l : iKratter et al.l l2010al ). Smaller instabilities may form brown dwarfs or giant 
planets f ^5 31). At late tirn es, several processes - including the clearing action of protostellar jets 
and winds ( Shu et al. 19871 ) - stop infall. Without a source of new material, the disk mass gradually 
declines with time. 

In general, all of the physical variables characterizing the cloud and the disk change with radius 
and time. However, Mj, M, and U often change slowly enough that it is useful to construct steady 
disk models with a constant M throughout the disk. Here, we develop the basic equations governing 
the evolution of the disk and use steady disks to show the general features of all circumstellar disks. 

An evolving gaseous disk sets the physical conditions in which small particles grow into planets. 
Physical conditions within the disk limit how planetesimals can form (21) and how solid planets grow 
out of planetesimals ( §5.1.5p . Once solid planets form, the gaseous disk provides the mass reservoir 
for giant planet atmospheres ( §5.2p and drives planet migration (see chapter by Morbidelli). 
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3.1. Basic Disk Dynamics 

To introduce basic concepts in disk dynamics we describe orbital motion in a gas disk and the 
radial flow induced by viscosity. The orbital velocity of the gas f<^, is set by the balance of radial 
forces - centrifugal, pressure and gravitational - as 

where p is the gas density and P is the gas pressure. Away from the immediate vicinity of pro- 
toplanets, gravity from the central star typically dominates. Even a self-gravitating thin disk has 
-^disk For P = Pq{R/Rq)~^ (the index n need only be locally valid), the orbital motion is 

/ nc2 \ 

v^ = VKil f , (3-3) 

in terms of the Kepler velocity, vk = \/GM^JR, and the sound speed, Cg = ^J^Pj with 7 the 
adiabatic index. Plausible disk models have an outwardly decreasing pressure {n > 0) almost 
everywhere. Thus pressure support typically gives sub-Keplerian rotation, but the correction is 
quite small, — vk ~ —Wi~^VK- It is often safe to ignore the pressure correction to orbital 



moti ons, but not when studying the motion of solids relative to gas (|Weidenschillinej|l977al : lYoudin 



20101, gm]). 

Accretion disks also have radial inflow, which is constrained by the laws of mass, angular 
momentum and energy conservation. Indeed accretion disks can be considered as machines that 
radiate energy as they transport mass inwards and angular momentum outwards. The viscous disk 
model offers the simplest means to understand how a disk manages this feat. Consider a thin ring 
with two adjacent annuli at distances R\ and R2 from a star with mass (Fig. [1]). Material orbits 
the star with angular velocities, Oi = \J GM^,/ Flf and = \/ GM^/i?2- If gas has viscosity, 
the differential rotation, Qi — ^2 < for R2 > Ri, produces a frictional (shear) force that attempts 
to equalize the two angular velocities. The resulting torques produce an outward flow of angular 
momentum. By moving a small amount of disk material onto distant, high-angular momentum 
orbits, large amount of mass can fall inwards to low angular momentum orbits. Mass accretion is 
biased towards inflow because specific (i.e. per unit mass) angular momentum can increase to very 
large values, but cannot fall below zero. The heat generated by viscous dissipation affects the disk 
temperature and the predicted spectra as described in ^3.31 

For a disk with surface density X! and viscosity mass c ontinuity and conservation of angular 



momentum lead to a non-linear diffusion equation for U (e.g.. iLvnden-Bell &: Pringldll974l : iPringle 
[iMll), 

The first term on the right hand side corresponds to viscous evolution; the second is a source 
term which is positive for infall from the cloud (i^ext = ^i) and negative for mass loss due to 
photoevaporation, disk winds, or accretion onto giant planets. Consider a simple model with = 



^The radial speeds associated with accretion produce negligible advective accelerations, Dvu/Dt ^ v\/R. 
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constant, Ucxt = 0, and X!{R,t = 0) = m5{R — Rq) /2ttRq. i.e. the initial mass m is in a narrow 
ring at radius Rq. The time evolution of the surface density is 



ni 



.-1 



(3-5) 



where the scaled distance and time are x = R/Rq and t = t/to = 12ut/RQ, and I1/4 is the modified 
Bessel function. Fig. [2] shows how this viscous solution asymptotically transports all of the mass to 
R = and all of the angular momentum to infinity. This evolution occurs on the viscous timescale 
r. 

We now consider the evolution of disks with a more general visco sity law, u = vniR/Rn)^ . Wit h 



a constant powerlaw /3, exact similarity solutions to eq. (j3-4p exist (|Lvnden-Bell &: Pringldll974l ). 
To develop intuition, we instead physically derive the approximate solution. After several viscous 
timescales have passed, the disk "forgets" the initial conditions (as seen in Fig. [2]). Conservation 
of angular momentum alone (without applying mass conservation needed to derive eq. |3-4j ) gives 



'dt 



1 



d 



i?3/2 QR 



MVR 
2tt 



'-uuVr 



(3-6) 



The two terms in square brackets represent the advection of angular momentum and viscous torques. 
For a steady disk with dU /dt = 0, the term in square brackets must equal a constant (independent 
of R). This constant represents the torque at the inner boundary, which can be neglected far from 
that boundary (as we show explicitly in eq. |3-17| ). Thus angular momentum conservation gives 



M 



oc R- 



(3-7) 



where the final proportionality assumes the powerlaw viscosity and constant M. From eq. (j3-ip 
the accretion speed follows as 

(3-8) 



VR 



3u 



To derive the evolution of disk mass (and M) in this limit, we adopt Rq as the outer edge of 
the disk. The outer radius changes on the local viscous timescale 



Rn 



ut 



1/(2-/3) 



(3-9) 



To conserve angular momentum and energy, the disk expands, requiring /3 < 2. Conservation of 
angular momentum, J ~ Mdisk^(^o)^0' drives the evolution of the disk mass as 



Afdisk 



The accretion rate through the disk is 



J 



VGM, 



-1/(4-2/3) 



^ ^ ^Mdisk ^ ^-(5-2/3)/(4-2/3) 

dt 



(3-10) 



(3-11) 



This physical derivation emphasizes the role of conservation laws in setting the global viscous 
evolution of the disk. The results are consistent with similarity solutions derived using Green's 
functions ( Lynden-Bell Sz Pringle 19741 ). 
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3.2. Transport Mechanisms and the a Disk Model 



The source of a physical mechanism to drive disk accretion is a longstanding problem. The 
molecular viscosity, z^moi = CgA — the product of sound speed, Cs, and collisional mean free path, 
A — is far too small. For the MMSN, A = ^imHCs/^^(JH2 ^ 0.8R^^^^ cm, where ^ = 2.4 is the 
mean molecular weight, niH is the mass of a hydrogen atom, and cj/Zj is the collision cross-section 
for the dominant constituent of the gas, H2. The resulting accretion timescale 



^acc,mol 



disk 



M 



i^mol 



7 X 10^^ 



R 
AU 



-4/7 



yr 



(3-12) 



vastly exceeds the age of the universe. Viable mechanisms for angular momentum transport are 
sometimes identified as "anomalous" sources of viscosity. Because long range interactions are often 
important, the analogy between a transport mechanism and the local viscosity is inexact. For 
sufficiently small scale fl uctuations, however, even self-gravitating disks are well approximated by 
the viscous prescription ( Lodato Sz Rice 20041 ) . 

Anomalous viscosity is most often associated with turbulence in the disk. Because the size 
of turbulent eddies can greatly exceed A, the effective turbulent viscosity can vastly exceed the 
molecular one, even for slow, subsonic turbulence. However, turbulence by itself does not explain 
accretion, velocity fluctuations must correlate for angular momentum to be transported outwards. 
Moreover it is difficult (if not impossible) for Keplerian shear to drive turbuler ice; the angular 
momentum gradient in disks is quite stable according to the Rayleigh criterion (IStone &: Balbus 
[l996). 

Over the past 50 years, theorists have considered a wide range of transport mechanisms: 
convective eddies, gravitational instabilities, internal shocks, magnetic stresses, orbiting planets, 
sound waves, spiral density waves, and tidal forces. Currently the "magnet o-rotational instability" 
(MRI) is the leading candid ate for a transport mechanism in low mass disks (jBalbus HawlevlEoEllI 



Papaloizou fc NelsonI lioo^l v In this mechanism, modest magnetic fields thread ionized material 



orbiting the central star. An outward (or inward) perturbation of fluid stretches and shears the 
magnetic fields. The resulting torque amplifies the original perturbation and, crucially, transports 
angular momentum outwards. Although this mechanism is attractive, the low ionization fraction 
of protostellar disks restricts the M RI to surface layers of the disk at many radii. Disks may 
then contain extensive "dead zones" (lGammielll996l ). where levels of transport and turbulence are 
reduced. In massive protostellar disks, g ravitational wa,yes an d gravitoturbulence are another likely 
source of angular momentum transport ( Lin fc Pringlel 1987 ). 

To sidestep fundai nental uncertainties of transport mechanisms, it is con venient to adopt a 
simple viscosity model ( Shakura Sz Sunyaev 19731 : Lynden-Bell Sz Pringle 19741 ). Setting u = aCgH 
- where H = Cgj O is the vertical scale height of the disk - leads to the popular "a-disk" model. The 
dimensionless parameter a < 1 (even <C 1) since large values would lead to rapid shock dissipation 
and/or gravitational fragmentation. Similar to the mixing length parameterization of convection, 
a-disk models allow progress despite ignorance of the underlying dynamics. Detailed simulations 
typically quote measured transport rates in terms of effective a values. 

This definition allows us to define t he three important timescales in a viscous disk (e.g., 
Lynden-Bell &: Pringlelll974l : iPring i3 ll98lh . The shortest disk timescale is the dynamical (orbital) 
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timescale, ~ ^ ^- 

The disk establishes hydrostatic equihbrium in the vertical direction on the same timescale ty ~ 

H/cs ^ n-^ ^ td. 

The disk cooling time, tc ^ U / D is the ratio of the thermal energy content (per unit area), 
U = CyST, to the energy generation or dissipation rate, D. For a viscous disk D = dvUQ'^/A 
(|Pringl(jll98ll ): C„ is the specific heat at constant volume. With C^T = (7 — 1) ^P/p for an ideal 
gas, the cooling time is 

In this classical result, the thermal timescale depends only on the local dynamical timescale, the 
dimensionless viscosity (a), and the equation of state (7). For a molecular gas, 7 ~ 7/5. The 
cooling time is a factor of roughly larger than the dynamical time: 

,,.„^„8,.„-. (JLy"\ ,3.15) 

The viscous timescale - t^, = E? /2>u - measures the rate matter diffuses through the disk. 
Using our expressions for the sound speed and the viscosity, t^ w {aQ)^^{R/H)'^ . Thus, the 
viscous timescale is 

Typically the disk is thin, H/R w 0.03-0.1; thus, the viscous timescale is 100-1000 times longer 
than the cooling time. 

The radial velocity in eq. (j3-8p becomes vr ~ acg H/R ^ 0.1a{H/Rfv^. With a < 1 and 
H/R < 1, is much smaller than both the orbital velocity and the sound speed. 



3.3. Viscously Heated Disks 

As shown above, t^ ^ tc > td ior a < 1, so the thermal properties of the disk adjust rapidly to 
changes in the surface density distribution. This property is very useful in describing the thermal 
properties of many astro physical disks , including active galactic nuclei, interacting binaries, and 
pre- main sequence stars ( Pr ingle 198ll ). We now describe the basic energetics of accretion disks 



with constant M 

Since gas remains on Keplerian orbits as it accretes, the specific energy release in falling from 
R + dR to R is GMi,dR/ i2R?)^ since half the potential energy goes into the increase in kinetic 
energy. The energy release per unit area for a disk accreting at M is Fk{R) = GMM^^/^iirR^). 
This result is completely independent of how the energy was released. 

The energy released by viscous dissipation does not simply match the local change in kinetic 
energy. To fully describe the energetics of steady viscous disks, we must keep the integration 
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constant when integrating eq. (|3-6p over radius to get, 

(3-17) 




The integration constant, Rj, represents the torque J exerted at the inner boundary, Rin as J = 
My/GMi,{y/Rin — \/ Rj) ■ The standard choice Rj = Rin is a zero torque boundary condition. 
Negative torques are not allowed for steady disks, as Rj > Rin would require I'U < 0. For Rj = 0, 
the maximum torque Jmax = M^/ GM^Rin matches the flow of angular momentum past the inner 
boundary. 



The laws of fluid dynamics in cylindrical coordinates (jShulll992l ) give the viscous dissipation 




dR) AttR^ 

As advertised, this expression does not simply match the local release of Keplerian orbital energy; 
far from the boundary {R » Rj), D{R) f» 3Fk{R)- Viscous disks transport energy (in addition to 
angular momentum) from the inner disk to the outer disk. Nevertheless, the rapid falloff with R 
means that most energy is dissipated close to the disk's inner edge. 

Now consider the total energy release from large R to the inner boundary. The Keplerian 
energy release is just Lk = GM^^M / {2Rin). The total viscous luminosity is 

L, = 2. r D(R)2.MR = { 1 - ''-M] (3-19) 

2 Rin y 3 V Rin J 

For the zero torque boundary condition [Rj = Rin) the luminosity simply matches the release of 
Keplerian energy. However the disk's luminosity increases due to work done by torques at the inner 
edge, up to Ld = "iLx for Rj = 0. For typical parameters in protostellar disks 

_ hGM.M ^ ( M \ fAU\ (R,n\-' 

"-'-^m-. o-i6/^^«|^io-B Mo yr-ij [wj fej ' 

where fd ranges from 1 (no torque) to 3 (maximum torque). 

The maximum disk luminosity occurs for a disk that extends to the stellar stellar surface, 
Rin = Ri,- The total accretion luminosity, 

_ f^GM^M 

^acc — ^ , (3-21) 

with 1 ^ /,t ^ 2, includes all the energy loss needed to come to rest on the rotating stellar surface. 
For a star rotating at breakup ~ 1. For a slowly rotating star the damping of the orbital 
kinetic energy gives twice the energy release ~ 2. Any difference Lace — > is emitted at 



Note that eq. 3.10 in IPringld (|l98ll l has a factor of two typo in the intermediate result (involving u) but reaches 
the correct final result (in terms of M). 
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the stehar surface. This difference must be positive (accretion should not cool the star), further 
constraining f^- As a consistency check, note that a disk with an inner boundary at the surface 
of a star, = i?^, that rotates at breakup, must satisfy the zero torque boundary condition to 
avoid Lace < Ld- 

In many cases, the ac creting star has a magnetosphere that truncates the disk at Rin > R-k 
(e.g. , iGhosh fc Lamblll979l ) . Material then flows onto the star along magnetic field lines, collimated 
onto hot spots, which are hot because the accretion energy is e mitted from a small fraction of the 
stella r photosphere. In most young stars, Rin 
20071 ) . Thus the hot spot luminosity 



3-5 Rj, (e.g., Kenvon et al. 19961 : Bouvier et al 



L 



hot 



Ld = La 



fdR* 
f*Rir 



(3-22) 



of the total accretion lumi nosity. For star that corotates with the disk's 



1 is expected (|Shu et al.111993 ;) 



can easily reach 60- 
inner edge, fd ~ 

To calculate the temperature structure of viscously heated disks, note that the upper and 
lower halves of the disk each radiate half of D{R). If the vertical optical depth r > 1, the disk 
photosphere then has effective temperature Tg, 



3GM^M 
87ri?3 




(3-23) 



Though Tf. is calculated as if the disk radiates as a blackbody, the disk's atmosphere will radiate as a 
stellar atmosphere with spectral lines, especially when > 1000-1500 K. The effective temperature 
declines as Tp cx i?"^/^ for R> a few i?^,: 



85 K 



M 



1/4 



10-« Mq yr- 



1/4 



-) 

Auy 



-3/4 



(1 - VRj/R) 



(3-24) 



In a simple, grey-atmosphere approach, the midplane temperature is a factor of r^/^ larger 
than Tg and is used to derive the scale height H, the viscosity z^, and other physical variables. 
More rigorous approaches calculate Tg and Td using a self-consistent prescription for the opacity 
throughout the disk. 

The temperature distribution in eq. (j3-23p allows us to derive the surface density of an a 
disk. Assuming the midplane temperature scales like the effective temperature (T^ ~ Tg oc R~^^^), 



v oc CgQ, oc r 



for the a prescription. Then eq. p-7p gives 



m) = To ( A 



-3/4 



(3-25) 



where Uq, the surface density at 1 AU, depends on the mass accretion rat e. Solving for the midplane 
temperature with realistic opacities yields E oc R-^ with j3 ^ 0.6-1 (e.g.. IStepinskill 19981 : IchambersI 
20091 ). 
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3.4. Steady Irradiated Disks 



Although viscous dissipation drives the overall evolution, radiation from th e central star also 
heats the disk ( Friedjung 19851 : Adams &: Shu 19861 : Kenyon Sz Hartmann 1987 ). If the disk is an 
infinite, but very thin, sheet, it absorbs roughly 25% of the light radiated by the central star. 
For a 1 Lq central star and disk accretion rates M < 10~^ Mq yr~^, emission from irradiation 
exceeds emission from viscous dissipation (eq. [3-20| ). If the disk re-radiates this energy at the 
local blackbody temperature, the radial temperature gradient of a fiat, irrad iated disk follows the 
gradient for a viscous disk, oc (|Friediung||l985l : lAdams fc Shulll986l l. 



D isks with vertical scale heigh t H absorb and re-radiate even more starlight ([Kenyon &: Hartmann 



19871 ) . IChiang &: GoldreichI (|l997l ) derive a general formalism for H and in a "passive" disk with 
negligible M. Defining 6 as the grazing angle that starlight hits the disk, the temperature of a disk 
that emits as a blackbody is 

1/4 

R 



(3-26) 



where is the stellar temperature. The grazing angle is 

Ri. d fh 



(3-27) 



where h is the height of the photosphere above the disk midplane. For a blackbody disk in vertical 
hydrostatic equilibrium, the grazing angle is the sum of a nearly flat component close to the star 
and a flared component far from the star: 



e ^ 0.005 



(R_ 



1 / R Y^' 

+ 0.05 — 



The disk temperature beyond a few tenths of an AU is then 

R \ "^/'^ / R^\ / T, 



(3-28) 



155 K 



AU 



Rq 



To 



(3-29) 



For M « 10 ^ Mq yr ^, the irradiated disk is roughly twice as hot as a viscous accretion disk. 

This temperature relation leads to a steeper surface density gradient in a disks. With u oc 
c^^^^ and cl oc T^, oc r^^/^^. Using this viscosity in eq. (l3-7p . 



-15/14 



where 



Si 



QArr 




(3-30) 



(3-31) 



The surface density gradient for an irradiated disk is steeper than the gradient for a viscous disk 
and is reasonably close to the gradient for the MMSN. 
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For identical a, hotter irradiated disks have larger viscosity and smaller surface density than 
cooler viscous disks. Integrating eq. (j3-3U|) over radius, the mass of an irradiated disk is 



When a ~ 10"^ — 10~^, this estimate is close to the observed masses of protostellar disks around 
young stars. 



3.5. Time Dependence 



Deriving more robust estimates of disk evolution requires a direct solution of eq. (|3-4|) . This 
exercise requires a model for the viscosity and a prescription for the thermodynamics and opacity 
of disk material. Analytic approaches assume a constant mass accretion rate through the disk. If a 
and r are simple functions of the local variables E a nd T, then the d i ffusion equation can be solved 



exactly for X'(t), M{t), and other disk properties (jStepinskil |199^ : Chamber^ 20091 ). Numerical 
approaches allow a and M to vary with ra dius. Some so lutions consider iterative solutions to 



the temperature and vertical structure (e.g. iHueso Sz GuillotI l2005l ) : others solv e for the vertical 



structure directly us ing techniques developed for the atmospheres of stars (e.g. 
D'Alessio etlllll998l ). 



BeU fc Lin 1994 



To compare these approaches, we consider a simple model for a viscous disk irradiated by a 
central star. We assume that the optical depth of cool disk material is dominated by dust grains 
with a constant opacity kq; warmer dust grains evaporate and have a smaller opacity: 



Kq 



(3-33) 



For material with roughly sola r metallicity, typical values are kq 

3. g., [chamber sl |2009l). With thi s opacity, we derive a self-consistent disk temperature and 



2 cmVg, T^^ap ~ 1380 K, and 



n = -14 (e.g.,IC 

scale height (as in Hueso &: GuillotI 2005) and solve for the t ime evolution of i7 using an explicit 
solution to the diffusion equation (as in Bath &: Pringk 19821 ). 

Fig. [3] compares analytic and numerical results for a disk with a = 10~^, initial mass M^ g 
= 0.04 M0, and initial radius Rq = 10 AU surrounding a star with M^^ = 1 Mq. The numerical 
solution tracks the analytic model well. At early times, the surface density declines steeply in the 
inner disk oc where dust grains evaporate) and more slowly in the outer disk {E oc R~^'^; 

where viscous transport dominates). At late times, irradiation dominates the energy budget; the 
surface density then falls more steeply with radius, S oc R^^ . 

Other approaches lead to similar time evolution in the surface density. Early on, a massive 
disk is dominated by viscous heating. For these conditions, the simple analytic estimate of the 
surface density yields E oc with n = 3/4 (eq. |3-25j ). close to results fo r the numerical 

solution (n = 0.6) and o t her analytic arid numerica.1 (n = 0.6 — 1) approaches (IBath &: Pringle 



19821 : iLin fc Pringi3ll99nl : Istepinskil[l998l : IChambersI [20091 : [Alexander fc Armitagell2009l ). As the 
disk ages, it evolves from a viscous-dominated to an irradiation-dominated system. Thus, the 
exponent n in the surface density relation approaches the limit (n = 15/14) derived in eq. (l3-30p . 
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Fig. m compares the evolution of the disk mass and accretion rate at the inner edge of the disk. 
In both solutions, the disk mass declines by a factor of roughly two in 0.1 Myr, a factor of roughly 
four in 1 Myr, and a factor of roughly ten in 10 Myr. Over the same period, the mass accretion 
rate onto the central star declines by roughly four orders of magnitude. 



3.6. Disk Instabilities and Fragmentation 



In addition to evolution on the viscous timescale shown in Fig. [3l all disks vary their energy 
output on much shorter timescales. In compact binary systems, these fluctuations range from small, 
10%-20%, amplitude flickering on the local dynamical time scale to large-scale er uptions, factors of 
10-100, that can last for several times the local viscous time scale (IWarnerlll995l). Although many 
pre-main sequence stars also display distinctive brightness variations ( Joy 19451 : Herbig 19621 ). the 
FU Ori variables provide the cleane st evidence for large-scale variations of the disk, rather than 
the environment or the central star ( Hartmann Sz Kenvon 19961 ). 

Theory suggests several types of instabilities in viscous disks (see Pringlel[l98ll ). In standard 
derivations of the structure of steady disks, radiative cooling balances heating from viscous stresses. 
However, radiative losses are set by local disk parameters; local parameters and an input accretion 
rate set viscous energy generation. Usually radiative losses can keep up with changes in disk 
structure; sometimes, radiation cannot balance viscous stresses, leading to a thermal instability. A 
limit cycle arises, where regions of the disk alternate between states where radiative losses exceed 
(and then fall below) the viscous energy input. This mechanism may produ ce FU Ori and other 
eruptions in the disks of pre-main sequence stars ( Hartmann fc Kenyon 19961 ). 

Viscous instabilities occur when changes in the local surface density do not produce parallel 
increases in the local mass transfer rate. From eq. (l3-7p . M oc I'E. In a steady disk, is fairly 
independent of U; thus, M changes in step with U. For the MRI viscosity mechanism, however, 
larger U leads to larger optical depths, less ionization, and smaller a. Thus, an MRI disk with 
growing (falling) surface density can produce a smaller (larger) viscosity, leading to an ever greater 
over- or under-density in the surface density. 

Although thermal and viscous instabilities change the temperature and surface density through- 
out the disk, they evolve on timescales much longer than the local orbital period. Massive disks can 
evolve on shorter timescales. If the local gravity in a region with size A overcomes rotati onal sup- 
port ( GE > n'^X) and thermal support (GL >.cl/X ] , this region can (begin to) collapse (|Safronovl 
1960l : iToomrel Il964l : IColdreich fc Lvnden-Belllll965l : IPaczvnskil Il978l ). Together these conditions 
require c^/{GE) < A < GS/Q"^. Collapse at any wavelength requires the disk satisfy the "Toomre 
instability criterion," 

LJ = 



Setting the disk mass Md ~ UR^ 



ttGE 
a stable disk has 



> 



< 1 



(3-34) 



(3-35) 



When the disk first forms, Md ~ iW*. Such "disks" cannot be thin because H/R ~ Cs/u^ > 1. 

I n a viscous ac cretion disk, the stability criterion can be re-written in terms of the accretion 
rate (jGammidlioOll ). With M = SttuU and i' = ac^Cl ^ an unstable disk has Mq 



> Sac^/G. To 



- 18 - 



From Disks to Planets 



Youdin & Kenyon 



evaluate the temperatures of unstable disks, we use Cg = {■jkT/ finiHy^'^ and set 7 = 7/5 and fi = 
2.4 for molecular gas: 

Mq > 2.4 10-^ qT3/2 Mq yr-i , (3-36) 

with T in Kelvins. For the observed accretion rates in very young stars, M ~ 10"'' Mq yr"-*^, 
unstable disks have aT^/^ < 0.4. If a is large (10~^), only very cold disks are unstable (T ~ 
10-15 K); smaller a (e.g., 10"^) allows instability in warmer disks (T ~ 50-60 K). 



4. FROM DUST TO PLANETESIMALS 

The accumulation of dust grains into planetesimals — solids greater than a kilometer in size — 
is the first step in the formation of terrestrial planets and giant planet cores. Several observational 
and theoretical reasons suggest the formation of planetesimals is a separate step. Observationally, 
remnant planetesimals in the Solar System and in extrasolar debris disks shows that growth some- 
times stalls before planets accumulate all planetesimals. Comets from the Oort cloud also suggest 
an intermediate stage between dust grains and planets (see Chapter by Moro-Martin) . 

Theoretically, the physical processes responsible for the growth of planetesimals differ from 
those relevant to the final stages of planet formation. As Section 5 describes, few-body gravitational 
encounters — both scattering and gravitationally focused collisions — establish the rates of growth 
for icy and terrestrial planets. By contrast, the sticking of dust grains involves electrostatic forces. 
During planetesimal formation, particle dynamics is qualitatively different. Drag forces exerted 
by the gas disk dominate the motions of small solids. Though not negligible, the drag exerted on 
km-sized or larger planetesimals is weaker than gravitational interactions (§5.2). 

While planetesimal formation is a common occurrence in circumstellar disks, understanding 
how it happens has proved elusive. Observations of planetesimal formation in action are indi- 
rect. Particles beyond cm-sizes contribute negligibly to images and spectra of circumstellar disks. 
Primitive meteorites record the conditions during planetesimal formation, but the implications for 
formation mechanisms are difficult to interpret — we need a better instruction manual. 

Especially beyond millimeter sizes, e xperiments show that particle collisions often result in 
boun cing or breaking instead of sticking ( Blum Sz Wurm 2000 : Zsom et al. 20 id : Weidling et al 



2OI2I ). Inefficient growth by coagulation is further complicated by the rapid infall of centimeter to 
meter sized solids into the star. These difficulties — often termed the "meter-sized barrier" — are 
explained in more detail in ^4.11 

Gravitational collapse is one way to overcome the growth barrier. The mutual gravitational 
attraction of a collection of small solids could lead to a runaway collapse into planetesimals — even 
when sticking is inefficient and radial drift is fast. While appealing, this path encounters theoretical 
difficulties when stirring by turbulent gas is included. Section 14.21 describes the current status of 
the gravitational collapse hypothesis. 

Even when self-gravity is weak, aerodynamic effects can concentrate solids in the disk. Particles 
tend to seek high pressure regions in the disk. This tendency causes the inward drift mentioned 
above. Particles can also concentrate in localized pressure maxima. Predicting the sizes and 
lifetimes of pressure maxima is a difficult (and currently unsolved) problem of disk meteorology. 

In addition to the passive response of solids to the gas disk, active particle concentration occurs 
when particles cause their own clumping by altering the flow of gas. Instabilities caused by gas 
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drag, notably the streaming instability, provide a clumping mechanism that is both powerful and 
amenable to study by direct numerical simulations. The strong clumping driven by the streaming 
instability is capable of triggering gravitational collapse into ~ 100 km planetesimals. Both passive 
and active particle concentration mechanisms are reviewed in ^4.3[ 

Theories based on complex non-linear dynamics must be tested against, and refined by, obser- 
vations. We discuss observational consequences of planetesimal formation models in §4.41 Unless 
stated otherwise, the numerical estimates in the section use the passively heat e d MM SN di sk (^3). 



see 



Chiang YoudinI (|20ld l and IVoudin 



For m ore detailed reviews of planetesimal formation, ^ ^ 

. j^or a thorough r eview of collision experiments and their relevance to planetesimal forma- 
tion, see lBlum fc WurmI (|2008l ^. 



4.1. The "Meter-Size" Barrier 

We discuss in more detail the two components of the "meter-size" barrier to planetesimal 
formation. The review of radial drift timescales in §4.1.11 also serves as an introduction to the 
dynamics of solids in a gas disk. The discussion of collisional growth and destruction in ^4.1.21 
couples dynamical models of collision speeds to the complex physics of contact mechanics. 



4-. 1.1. Radial Drift and the Basics of Disk Aerodynamics 

The aerodynamic migration of small solids imposes the most stringent timescale constraint 
on planet formation: ~ 100 years. Aerodynamic radial drift arises because solids encounter a 
headwind as they orbit through the gas disk (eq. |3-3| ) . This headwind removes angular momentum 
from particle orbits, causing their inspiral. Infall speeds are fastest for solids near roughly meter 
sizes. The critical size is actually below a meter in standard disk models — especially in their 
outer regions. So the "meter-size" barrier is a slight misnomer, but it has a better ring than the 
"millimeter-to-tens-of-centimeters-size" barrier. 

Radial pressure gradients in gas disks set the speed of the headwind. Plausible disk models are 
hotter and denser in the inner regions; on average, the radial pressure gradient is directed outwards. 
If the radial pressure gradient is directed inwards, a tailwind — and outward particle migration — 
results. 



We express the headwind speed as the difference between the Keplerian velocity, fx = ^/GM^JR 
f2R, and the orbital speed of the gas, Vg^^: 

dP/dlnR RV'^^ _i 

r?UK = f K - Vg^^ ^ 25 — - m s , (4-1) 



2pgVK V^U^ 

where P and pg are the pressure and density of the gas, and ?7 ~ Cs/f|; ~ [H/R)'^ ~ 10~^ is a 
dimensionless measure of pressure support. In disks hotter than our passive model, headwinds and 
drift speeds are faster. 

To derive eq. (j4-ip . compute radial force balance assuming (correctly) that the radial pressure 
acceleration, fp^R = —p~^dP/dR, is weak compared to the centrifugal acceleration. Equivalently 
we can reproduce eq. (|4-ip by balancing the pressure and Coriolis forces, fp^R + fcor,R = with 
/cor,Ji = -2QriVK- 
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Drag forces set the response of particle orbits to the gas headwind. We express the drag 
acceleration felt by a solid particle as 

/d,ag = -A^/ts , (4-2) 

where Av is the particle velocity relative to the gas and ts is the aerodynamic damping timescale 
for this relative particle motion. 

The value of tg depends on particle properties - such as the internal density, p,, and spherical 
radius, s — and on properties of the gas disk — pg and Cg - as: 

tfP = p.s/{pgCs) if s < 9A/4 (4-3a) 

tf°^^' = tfP ■As/{9X) if 9A/4 < s < A/(4Ma) (4-3b) 

• (s/^f^^ if A/(4Ma) < s < 200 A/Ma (4-3c) 

4"'^ = • 6/Ma if s> 200 A/Ma (4-3d) 

where Ma = |At;|/cs, A oc l/pg is the gas mean free path, and Re = 4sMa/A is the Reynolds number 
of the flow around the particle. The cases are written in order of increasing particle size: Epstein's 
Law of drag from molecular collisions. Stokes' Law for viscous drag when Re < 1, an approximate 
intermediate Re case, and the drag from a fully developed turbulent wake for Re > 800. The 
turbulent drag force is more relevant for fully formed planetesimals and is commonly expressed as 

^ \Av\ Cd 2 I a i2 / a .\ 

i^drag = -m-^ = ^TTS Pg|A^;| , (4-4) 



where the drag coefficient, Cd ~ 0.44 (|Adachi et al.lll97()l : FvV^idenschillinglll977al V 



The dynamical significance of drag forces is measured by comparing the stopping time and the 
orbital frequency, via the parameter 

Ts = nts . (4-5) 

For Tg ^ 1, particles are carried along with the gas; for Ts S> 1, gas drag is a small correction to 
Keplerian orbits. Fig. [5] plots Tg for a range of particle sizes in our passively heated MMSN. At 
least in the inner disk, objects near meter-sizes have Tg ~ 1. In the outer disk, where gas densities 
are lower, smaller solids have the critical Tg = 1. As we now show, Tg = 1 solids have the fastest 
radial drift speeds. 

To derive the particle drift caused by the gas headwind, we consider the equations of motion 
for a particle in cylindrical coordinates, r and (j), 

R-R^"^ = -vl/R - R/tg (4-6) 
R4> + 2R^ = -{R^ - Wg,</,)/ts • (4-7) 

To find the steady drift solutions, we make several approximations that are valid when drag forces 
are strong. We neglect the radial inertial acceleration, i?, and express the azimuthal motion as a 
small deviation from the Keplerian frequency, (j) = Q + 5v^/ R where \5v^\ <^ QR. The azimuthal 
acceleration is then cj) k, Q k. —2>f2R/ {2R). 

The radial drift speed follows from equations ()4-ip . (I4-5|) . (I4-6|) and ()4-7p as 
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Solids with Tg = 1 have the fastest infall speed, —R = r]VK- The corresponding timescale 

min(t/j) ~ iv^y^ ~ 200(i?/AU)^3/^^ yr (4-9) 

is a very strong constraint on growth. This constraint is the main element of the meter-sized growth 
barrier. Fig. [6] plots the radial drift timescales for a range of particle sizes. 

To complete this brief introduction to particle aerodynamics, we give the azimuthal drift speed 
of solids through the (sub-Keplerian) gas as 

R4> - Vg,(j, = SVfp + r]VK = rjVK i _J 2 ' (4-10) 

Large, Tg ^» 1 solids experience the full r/ux headwind, yet their radial drift is slow because their 
inertia is so large. Small, Tg <C 1 solids are dragged by the gas and only feel a mild headwind. We 
thus see why radial drift is fastest near Tg « 1. For these intermediate sizes, drag forces are strong 
enough to overcome particle inertia, but not so strong as to cause perfect coupling. 

These idealized calculations explain the basics of radial drift. A pressing question is whether 
ignored effects could mitigate the radial drift problem. The existence of a headwind is the most 
crucial assumption, and it can vanish in localized pressure maxima as addressed below. Even 
when these maxima exist, headwinds still prevail in the majority of the disk. We also assume 
that aerodynamic drag only affects the solids, and not the gas component of the disk. When 
the distributed mass density of solids pp becomes comparable to t he gas density Pp - , then it is no 
longer acceptable to ignore the feedback of drag forces on the gas. Nakagawa et al. ( 19861 ) showed 



how drift speeds become slower when /Op > Pg. This feedback is also the source of powerful drag 
instabilities — both shearing and streaming — that we address below. Thus there is no simple way 
to ignore the radial drift problem — its resolution has consequences for how planetesimals form. 



4-1.2. Early Collisional Growth 
Planetesimal formation begins with the collisional agglomeration of dust grains into larger 



solids, a process that is observed to proceed up to mm-sizes in T Tauri disks (IWilliams &: Cieza 



20111). The conceptually simplest mechanism to form planetesimals is for this collisional growth to 
proceed past kilometer sizes. However both direct experiment and theoretical arguments show that 
coagulation beyond mm-sizes is inefficient at best. This inefficiency is particularly problematic due 
to the timescale constraints imposed by radial drift. The combination of inefficient sticking and 
rapid infall together comprise the formidable "meter-size" barrier. 

Although collision rates do not rule out rapid growth, they place tight constraints on the 
sticking efficiency To make this conclusion, we approximate the mean collision time as tcoii ~ 
l/(i7r), where r ~ Sp/{p,s) is the vertical optical depth. This approximation for the collision rate 
is good when Tg ^ 1, and we show below that it suffices for Tg = 1. The ratio of collision to drift 
timescales for Tg = 1 solids is thus roughly 



idrift 



r? f R \^/V0.01', 
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where we assume Epstein drag, appropriate for the outer disk. We use the result of hydrostatic 
balance that ~ p^c^jQ. When the collision time exceeds the drift time, collisional growth is 
ruled out. Even when the two are close, growth requires an efficient rate of sticking per collision. 
This constraint is most severe in the outer disk. 

While turbulent motions increase collision speeds, they do not increase the collision rate above 
the geometric estimate used to derive eq. (14-lip . The reason is that turbulence also increases the 
particle layer thickness, i?p, thereby decreasing the mean particle density, pp. We can compute the 
collision rate due to turbulence as t~^y^ ~ nav. Here the particle number density is n ~ ^/{Hpirip 
where rup is the particle mass. The particle layer thickness due to turbulent stirring i: 



Hp — Ha 




(4-12) 



This well-known result (jCuzzi et al.l Il993l : ICarballido et al.l bood : lYoudin k Lithwickl hopj ) nor- 
malizes the turbulent diffusion, D, to the dimensionless parameter, od = D/{csHg). The cross 
^ and the relative velocity du e to turbulent motions is, v ~ -^Z UTg/ (1 + t^)cs 



section is a 



dMarkiewicz et al.lll99ll : IChiang k Youdinlboid ). For Ts > 1, the collision rate necessarily agrees 
with the optical depth estimate. For < 1, the collision rate is also independent of the strength of 
turbulence as t~^^^ ~ Zfi. These cases agree at Ts ~ 1 and confirm the constraint set by eq. (j4-lip . 

Collision rates are not the only concern. Collisions can also result in bouncing or fragmen- 
tation that stalls, or even reverses, growth. Below speeds of ~ 1 m s~^, small dust grains stick 
efficie ntly as a consequence of van der Wa,als iri teractions and the efficient dissipation of kinetic en- 
ergy ( Chokshi et al. 19931 : Blum fc Wurm 200d ) . As particles grow and as collision speeds increase. 



the collisional kinetic energy increases. Short rang e sticking forc es cannot match this increase in ki 



"J — ^¥ — I 

netic energy, because they are surface area limited (lYoudinll2004l ) . Experimental wo rk confirms that 



mm-sizes (Blum fc Wurm 



collis ions between equal mass objects do not produce growth beyond 
20081 ). 

Collisions between lower mass projectiles and higher mass targets offer another route to growth. 
In this scenario impact speeds exceed the meters-per-second value expected to produce growth. As 
explained above, since small solids are tied to the gas flow, they impact larger solids (which decouple 
from the gas) at the full headwind speed, rjVK ^ 25 m s~^. Indeed when the latest experimental 
results are combined with dynamical estimates of collis ion speeds for a dispersion of particle sizes, 
growth stalls at only millimeter sizes (|Zsom et al.ll20ld ). 

Based on obs erved SEDs, disks likely find a way to overcome these obstacles and grow solids 
beyond mm-sizes (jWihiams fc Ciezall201lh . The mechanisms responsible for enhanced coagulation 
remain unclear. The particle con centration mechanisms discussed below could augment particle 
sticking. As shown explicitly in I.Tohansen etHI (l2009bl ). collision speeds are reduced in dense 
particle clumps. 

Most experimental work on grain-grain collisions uses porous silicates. If ices are stickier, 
growth beyond mm-sizes is possible. In the low pressure of disks, ices sublimate; there is no liquid 
available to make the equivalent of wet snow. Saturn's rings are an excellent labo ratory to explore 



the outcomes of gentle, ~ mm s ^, collisions between ices (lYoudin k Shull2002l ). Here, sticking 



'To accommodate the even mixing of small grains (not our current concern) we require Hp < Hg 
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forces are constrained by their inabihty to overcome tidal shear and produce growth beyond ~ 5 m 
objects. Terrestrial experiments on low temperature ices sug gest that cm-sized frosty objects stick 
at collision speeds below ~ 0.1 m s~^ ( Supulver et al.l 119971 ). While possibly a crucial ingredient, 
this limit appears insufficient to allow icy surfaces to bridge the meter-size barrier. 



4.2. Gravitational Collapse of Solids into Planetesimals 



Self-gravity provides a qualitatively different route to the formation of planete s imals . In- 
stead of bottom-up growth, the gravitational instability (GI) hypothesis of Safronov ( 19691 ) and 
Goldreich fc Ward! ([l973) offers a top-down approach. In this theory, a sea of small solids collapses 
coherently into a gravitationally bound planetesimal. This collapse does not rely on sticking forces, 
proceeds faster than radial drift, and bypasses the meter-size barrier. 

The gravitational collapse hypothesis encounters several theoretical di fficulties. The crucial 
issue is the ability of turbulent gas to prevent collapse ( Weidenschilling 19951 ). Until recently, these 
theoretic al obstacles seemed insurmountable. Progress in coupled particle-g a .s dynamics ha s led to 
a revival (I Youdin fc Shull2002l : Ijohansen et"alll2006l . l2007l : ICuzzi et ahlboosl : lYoudinllioilal ). Some 
of these mechanisms use aerodynamic concentration as the initial concentration mechanism ( 94. Sp 
but all eventually rely on self- gravity for the final collapse to solid densities. 

We focus in this subsection on "pure" gravitational collapse from a relatively smooth back- 
ground. Although separating these processes from aerodynamic concentration is artificial, this 
historical approach allows us to isolate the main issues of each mechanism. We first discuss the 
standard model of gravitational collapse of a disk of solids, which has many similarities to gravita- 
tional instabilities in a gas disk (§5). We then briefly describe how gas drag changes this standard 
picture, a research area where progress is still being made. 

The simplest criterion for gravitational collapse requires self-gravity to overcome the tidal 
distortion of the central star. This condition is met when the particle density exceeds the Roche 
limit, 



/Op > /OR - 0-6^ ^ 130- 



R 
AU 



-3/14 



i?3 F 

M^/Mq and F is the mass enhancement factor for the MMSN from §2. 



(4-13) 



Sekival (|l983l l 



where m- 

derives this result for the case of a disk midpl ane with solids perfectly coupled to an incompressible 
gas, making use of th e powerful form alism of lGoldreich k Lvnden-Beiil (|l965l ^. The relation of this 
cri tical density to the Toonirel ( 1964) Q criterion for GI is discussed (in the context of planetesimals) 
bv lChiang YoudinI (|20ld ) and lYoudh] (|2011al ). 



Mr 



ptml 



For a given particle surface density, eq. (|4-13p implies that planetesimals form with a mass 
^p/ Pr- After contraction to solid densities, the planetesimal size is 



-Rptml 



1/3 2/3 
P' PR 



AU 



km . 



(4-14) 



Though the current relevance is not so clear, this kind of estimate played a key role in defining the 
canonical planetesimal size to be near a kilometer. 

To satisfy the density criterion of eq. (j4-13p . solids must settle vertically to a midplane layer 
with thickness i?R = Up/pn. Even the faintest whiff of turbulence probably produces a much 
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thicker layer. A lthough the disk midplane could be a "dead zone" devoid of magnetized tur bulence 
(lGammielll996l l. interactions among particles can drive enough turbulence to halt settling ([Youdinl 
2O10l ). 

Vertical shear inst abilities usually prevent the sedimentation of small particles into a layer 
thinner than i?^ ~ r]R dWeidenschinindllQSOl : lYoudin fc Slmlbood ). As particle inertia in the mid- 
plane increases from sedimentation, solids begin to drag the midplane gas towards the full Keplerian 
speed. As in the Kelvin-Helmholz instability, the vertical shear with the overlying particle-poor 
gas drives overturning. With Hrj/H^ ~ 200, GI seems ruled out. 

The revival of the GI hypothesis requires abandoning two faulty assumptions. The surface 
density of solids can increase above MMSN — or any initial — values. The evolution of solid 
and gas components decouples due to drift motions (IStepinski fc: Valageas The radial drift 

of small solids from the outer disk generically leads to "particle pi l eups," a snowplow effec t that 
increases the surface density in the inner disk dVoudin fc Shulbooi lYbudin Chiang|l2004l ). The 
local concentration mechanisms discussed in §4.31 can be even more powerful. 

The critical Roche density is also too stringent. Planetesimal formation can be triggered 
when pp > pg, a criterion about a hundred times less severe than the Roche limit in eq. (14-13^ . 
Several interesting effects arise when the particle density approache s the gas density. Vertical 
shear instabilities lose their ability to overturn a layer that is so heavy ( Sekiya 19981 : Youdin Sz Shu 
2OO2I ). W hen disk rotation is included, the case for particle inertia halting vertical overturning is 



less clear ( Gomez Sz Ostriker 2005 : Lee et al. 2O10l ). However, when perfect coupling is relaxed, 
and streaming instabilities appear, the relevance of pp > pg reemerges as the threshold for strong 
clumping, as described below. 



When gas drag is present, the Roche density is not the relevant criteria for GI. WardI ( 19761 . 
2000l ^ investigated a dissipative mode of GI that has no formal stability threshold. Collaps e always 
occurs in principle, but it becomes slower and spans a larger radial extent for small particles. lYoudin 
(|2011al ) included radial turbulent diffusion, and showed that radial spreading — not vertical stirring 
— is the dominant stabilizing influence for dissipative GI. When vertical stirring is accounted for, 
it turns out that pp > pg is typically required for dissipative GI to proceed faster than radial 
drift. The result that dissipative GI depends so simply on particle inertia is m ostly a riumerica l 
coincidence and relies on the fact that Qg^/2rj ~ 1 in the MMSN, see eq. (55) of I Youdinl (l2011al ). 
The important point is that the relevance of particle inertia — specifically the pp > pg criterion — 
has been established for a range of mechanisms. 

Although this lesser degree of particle settling is substantial, it may require a local enrichment 
of the disk metallicity, Z = Ep/Ug. If the particle scale- height is set by particle-driven turbulence 
to H„, then the particle density exceeds the gas density if 



Z > 



r]R 



0.014- 



Auy 



2/7 



(4-15) 



again with m^, = M^/Mq. The near agreement with Solar abundances is remarkable and coul d 
be related to the correlation of giant planets with host star metallicity (lYoudin &: Shul |2002|) . 
Assu ming that the stellar photospheres reflect the abundance of solids in the disk (i Fischer &: Valenti 
2OO5I'). the early format ion of planetesimals could be a crucial factor in the formation of gas giants 
(jjohansen et al.ll2009bl l. 
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Though poorly constrained, the role of external (not particle-driven) midplane turbulence may 
be interesting. For small solids, constraints on the level of turbulence that allows settling to pp > pg 
are quite stringent. Using eq. (|4-12p . sedimentation to pp > pg requires that midplane turbulence 
satisfy 

aD<27rZ\,^W-^Zl,4-{^^] . (4-16) 



cm V 10 AU 



Thus when trying to form planetesimals via GI it helps to have some combination of weak turbu- 
lence, particle growth ar id enriched me tallicity (Z'p/Z'g). These requirements become more stringent 
towards the inner disk (|Youdinll2011al ). 

Thus even in the GI hypothesis, particle growth by coagulation plays a crucial role. Particles 
must grow until they decouple from the gas. Provided this growth occurs, GI — likely aided by 
other concentration mechanisms — provides a plausible way past the meter-size barrier. 



4.3. Aerodynamic Particle Concentration 



We now consider aerodynamic processes that can concentrate particles even when self-gravity is 
negligible. Many of these processes rely on the presence of turbulence in the disk. This connection 
raises a general question: does turbulence help or hinder planetesimal formation? By stirring 
particles, turbulence increases their collision speeds which can lead to more destructive collisions. 
Furthermore, the diffusive effects of turbulence oppose particle settling and concentration. On 
the other hand, turbulence can concentrate particles in a variety of ways. Which tendency wins 
depends on details, notably particle size. Since small solids with Tg ^ 1 drift and settle slowly, 
they require much weaker turbulence to participate in aerodynamic concentration. 

Localized pressure maxima are very powerful particl e traps. When t he pressure bump takes the 
form of an axisymmetric ring, the trap is very effective (|WhiDDlelll97i l. Solids migrate into these 
rings and accumulate at the center where they encounter no headwind. The MRI naturally produces 
axisymmetric pressure bumps, via the generation of zonal flows that are sorn ewhat analogous to 
the surface winds of Jupiter (jJohansen et al.l l2009al : iFromang &: Stonel bnOflh . The relevance of 
MRI-induced pressure maxima is subject to two caveats: turbulent stirring associated with MRI 
may lead t o destructive collisio ns and the disk midplane may be insufficiently ionized for the MRI 
to operate ( Turner et al. 20ld ). 

Non- axisymmetric pressure maxima can also trap particles. When the disk is young and 
mass ive, spiral arms in the gas probably provide an important source of turbulence (iRice et al. 
20061 ). However, this phase of disk evolution may be too turbulent and/or brief for significant 



plane tesimal formation. Isolated pressure maxima take the form of anticy clonic vortice s (jChavanis 
2003). Vortices are embedded in, and thus flow with, the sub-Keperian gas (jYoudinlboiol ). Although 
the vortex center is not a stable point for particle concentration, a p oint ups t ream (in the direction 
of orbital motion) is. The implications for vortex size i s discussed inlYoudin The formation 

and survival of vortices is a topic of ongoing research ( Lithwick 20091 ). 

We have so far focused on particle concentration over many orbits where disk rotation and 
Coriolis forces play a central role. Small turbulent eddies have short turnover times, teddy <C 1/i?, 
and are unaffected by rotation. In this regime pressure maxima occur not at the centers of anti- 
cyclonic vortices, but between vortices of either sign. The concentration of heavy particles in these 
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regions of low enstrophy (vorticity-squared) was first described in the fluid dynamics community 
(|Maxevlll987l). 

et al.1 (|200lh applied small-scale concentration to protoplanetary disks. They showed 
that ~ 1 mm solids — specifically the chondrules that are discussed in ^4.41 — can concentrate at 
the "inner" or dissipation scale of turbulence. These are the smallest eddies that have the shortest 
turnover time tj. Particles with a matching stopping time, ts)~ ~ 30 s, are preferentially flung 
from the eddies and concentrated between them. Chondrules can plausibly satisfy this condition. 
The ability to concentrate such small particles makes this mechanism unique. The relevance of 
such brief concentrati ons is unclear. The cha racteristic mass involved is also quite small, at most 
that of a 10 cm rock dChiang k Youdinlboid ). 

To overcome these issues, Cuzzi etliD ( 20081 ) developed a model that concentrates chondrules 
on larger scales that contain enough mass to form ~ 100 km-planetesimals. This model involves a 
somewhat speculative extrapolation. In particular, it assumes that all scales of a turbulent cascade 
contribute equally to the concentration of chondrule-sized particles. This assumption is a signifi- 
cant deviation frg r n the original mechanism that requires eddy and stopping times to match. See 
Chiang Youdin for further discussion, which concludes that more study of this intriguing 



mechanism is required. 

Clearly particle concentration mechanisms are fraught with uncertainties in the detailed dy- 
namical behavior of gas in protoplanetary disks. Some — certainly not all — of these uncer- 
tainties are overcome by the realizat ion that particles caii cau se their own concentration by col- 
lectively altering the gas dynamics ( Goodman Sz Pindor 2000l ). In the streaming instability of 



Youdin &: GoodmanI (j2005l ). particle concentrations arise spontaneously from radial drift motions. 



As described in §4.1.11 these drift motions are an inevitable consequence of pressure support in 
disks. The linear growth of streaming instabilities is strongest for pp > pg, because particle in- 
ertia must be large for drag feedback to ii ifiuence gas motions . Wh en pp < pg, growth is fastest 
for Tg ~ 1, when drift speeds are fastest ( Youdin Sz Johansen 200?! ). While streaming instabili- 
ties involve complex dynamics — 3D inotions of both the gas and solid components — simplified 
toy mod els (Goodman &: Pindor 2000 : Chiang &: Youdin 20 id ) and considerations of geostrophic 



balance (jJacquet et al.ll201ll ) help explain how particle density perturbations self-reinforce. 

Numerical simulations show that the non-linear clurnping from the streaining instability can 
be qui te strong (jjohansen &: YoudinI 120071 : 1 Johansen et al.ll2009bl : iBalsara et al.ll2009l : iBai &: Stone 



2010al V Particle densities > 10 pg are achieved in the absence of self-gravity, and clumping tends to 



increase with numerical resolution. The conditions for strong clumping are similar to those giving 
rapid linear growth: partial decoupling, Ts > 0.1, and large particle inertia pp > 0.2pg. 

When vertical stratification is included, the midplane particle density evolves consistently due 
to settling and stirring by both streaming and vertical shearing instabilities. In th ese simulation s, 



there is a critical disk meta llicity for particle clumping which is slightly super-Solar (jJohansen et al 



2009bl : iBai fc Stonell2010al ). consistent with eq. (j4-15p . This metallicity vari e s with the ra,dial pr es- 



sure gradient, rj; smaller gradients promote clumping (jJohansen et al.l 120071 : iBai Sz Stonell2010bl ). 

In §4.21 we noted that GI depends on particle growth by coagulation. Since particle sedimen- 
tation to Pp ^ Pg is a crucial prerequisite, particle growth remains essential when the streaming 
instability provides the initial particle concentration. However growth need not result in a single 
particle size, or a very narrow size distribution. Though the smallest solids participate less in 
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clumping by streaminK instabili ties, including a dispersion in particle sizes does not prevent strong 
clumping ( Johansen et al. 2007I ). 

The particle concentrati on produced by the st reaming instability is more than sufficient to 
trigger gravitational collapse. Johansen et al. ( 200?! ) formed gravitationally bound objects equiva- 
lent to ~ 500 km planetesimals within only a few orbits of initial collapse. More re cent simulations 
sugges t the formation of lower mass objects with equivalent ~ 100 — 200 km radii (jJohansen et al 



2009bl l. The crucial differences are the inclusion of the MRI in the earlier study and smaller particle 
sizes in the second. A more thorough investigation of parameter space, combined with resolution 
studies, is required. These results exceed the standard estimate of km-sized planetesimals because 
gravitational collapse occurs not from a smooth background, but from aerodynamically concen- 
trated clumps. 



4.4. Observational Constraints on Planetesimal Formation 



We now discuss how observations constrain dynamical theories of planetesimal formation. 
The Solar System provides the most detailed information on planetesimals, and allows comparison 
between the inner asteroidal reservoir and the Kuiper belt objects and comets of the outer Solar 
System. The crucial issue is to what extent today's planetesimals reveal the clues of their formation, 
especially after ~ 4 Gyr of dynamical, collisional and thermal evolution. 

Primitive, undifferentiated meteorites give us a hands-on view of the composition of planetes- 
imals. The most common of these are the aptly named "ordinary chondrites." With filling factors 
up to 90%, they are primarily composed of 0.1-1 mm ch ondrules. Chondrules are glassy spheres, 
poetically referred to as "fiery drops of rain" (ISorbvlllSed ^. The o rigin of c hondr ules — in particular 
their source of heating — is debated and beyond our scope, see iHewinj (|l99(il ). The prevalence of 
chondrules in ordinary chondrites s trongly motivates furthe r investigation of the mechanisms that 



could concentrate solids this small ( Cuzzi et al. 2001, 20081 ). 



Despite this attractive conclusion, chondrules may not be the universal building blocks of all 
planetesimals. Because their abundances mo st closely match Solar, the CI class of chondrites is 
considered the most primitive (Lodders 20031 ). Yet CI chondrites contain no chondrules. It is also 
unclear whether chondrules were present in the first generation of planetesimals. Most chondrules 
formed at least 1.5 Myr and up to 4 Myr after the rarer CAIs (Calcuim- Aluminum Inclusions; 
Connellv etaPboOTl : iKrot et al.l[2007l l. Thus, planetesimals probably formed before major chon- 
drule forming events, especially the planetesimals that formed the cores of Jupiter and Saturn. 
Since planetesimals that form early will trap more radioactive heat and differentiate, i t seems likely 



that the undifferentiated chondrites represent a later phase of planetesimal formation (jKleine et al 



20051 ). The relation between chondrules, meteorites and planetesimal formation continues to be the 
focus of intense interdisciplinary research. 

Planetesimals that remain in the asteroid belt can also provide clues to th eir formation. The 
radia lly banded zonation of different spectral classes of asteroids is well known ( Gradie Sz Tedescol 
1982 ). This observation suggests separate forma tion epochs, wi th each event creating a "clan" of 
chemically and spectrally similar planetesimals. Q ffl) proposed large-scale, drag medi- 
ated GI as the cause of these events. 

The size distribution of objects within planetesimal belts provides other clues to their for- 
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mation. Breaks in the size distribution — i.e. changes in its powerlaw slope — point to shifts 
in formation an d /or ero sion processes. The asteroid belt has a break near a radius ~ 50 km. 
Morbidelli etaP (j2009al l argue that the asteroids with radii > 50 km reflect their initial sizes. 



Specifically they assert that the largest asteroids have undergone minimal collisional evolution and 
could not have formed via collisional growth of smaller planetesimals. Since most of the mass is 
contained in the largest asteroids, their model plausibly produces the numerous small objects below 
the break via collisional disruption. That interpretation places GI as the preferred forma tion mech- 
anism. By including streaming instabilities, the simulations of iJohansen et al.l (120071 ) predicted 
that large initial sizes were possible. C onclusively proving tha t a size distribution is unobtainable 
by collisional growth is rather difficult. Weidenschilling ( 2O10l ) contends that collisional growth of 
asteroids can be accomplished starting with 0.1 km planetesimals — which themselves presumably 
grew by coagulation past the meter-size barrier. 

^Curiously, the Kuiper belt also has break in its size distribution at ~ 50 km radii (iBernstein et al. 



20041 ). This break is not measured directly; a combination of an observed lum inosity distribution 



and an estimate of the albedo yields the distribution of radii ( Petit et al. 20081 ). Ongoing surveys 
of the Kuiper belt seek to provide constraints on the size distribution for the various components 
of the Kuiper belt. 

Understanding the origin of the break requires a model for KBO formation and dynamical inter- 
actions with gas giants. Reproducing the observed size distribution with collisional growth models 
requires an initially massive Kuiper belt followed by dynamical depletion; a break occurs when 
deplet ion excites erosive collisions among KBOs with radii below the break ([Kenvon &: Bromlev 



2004d ) . The break radius depends on excitation; more (less) excitation by more (less) massive gas 



giants yields a break at larger (smaller) radii. Matching the location of the break and the apparent 
slope of the KBO size distribution below the bre ak requires numerica l calculations with grow th and 
depletion, which are an active area of research ( Kenyon et al. 2008a : Morbidelli et al. 20081 ). 

The r elative roles of coll isional and dynamical depletion affect the interpretation of the KBO 
size break. Pan &: Sari (j2005l ) argue that the break is not primordial but due to ongo ir ig col lisional 
erosion that continues to push the break to larger sizes. However iNesvornv et al.l (|201ll ) claim 
that this collisional history is ruled out on two grounds. First the collisional strengths required for 
such destruction are too weak. Second, such an intense collisional bombardment would destroy the 
observed Kuiper belt binaries. 

The observed binary fraction in the cold, classical Kuiper belt is > 20% ( Noll et al. 20081 ). The 
colors of the two components of Kuiper b elt binaries are nearly identical, a fact interpreted as repre- 
senting a common chemical composition ( Benecchi et al.ll2009l ). This observation provides the most 



compelling support for the GI hypothesis in the outer Solar System, or perhaps anywhere. Gravita- 
tional collapse can naturally produce binary planetesimals as a consequence of angula r momentum 
conservation during the contraction of a swarm of small solids ( Nesvorny et al. 2010l ). Binaries — 
and higher-order multiples — formed this way should have the same chemical composition since 
they formed from the same well- mixed clump of small solids. While mechanisms for the dynamical 
capture of KBO binaries are well developed ( Goldreich et al. 20021 : Noll et al. 20081 ). these models 
do not obviously explain matching colors. Moreover the ph ysical conditions requ ire for capture, 
make the colli sional survival of these b inaries questionable ( Nesvorny et al. 201 ll ). especially for 
wide binaries ( Parker &: Kavelaard 20 id ). 
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Outside the solar system, exoplanets and debris disks inform the prevalence and consequences 
of planetesimal formation. The higher incidence of giant planets around stars with super-Solar 
metallicities (discussed in ^2.3p might be tied to planetesimal formation. As shown in eq. (|4-15p . this 
connection is strongly suggested by the super-Solar disk metallicity threshold fo r strong clumping 
by st reaming instabilities. Since the disk metallicity can increase over time (jYoudin Chiang 
20041 1. this threshold does not imply that lower metallicity stars can never form planetesimals. 

Indeed, the streaming instab ility/GI model explains why lower metallicity and lower mass stars 
should form less massive planets ( Johansen et al. 2009bl ) . Either directly or by the passage of time, 
enriching the disk metallicity involves the loss of gas. Thus the initially low metallicity systems 
that require enrichment are less likely to form giant planets. This conclusion is especially true in 
the lower mass disks t hought to surround lower mass stars. These general trends are revealed by 
radial velocity surveys dSousa et al.ll2008l : Ijohnson et aDl20ld l . The Kepler tr a,nsit survey will test 



these trends, s ince it is finding striking numbers of small, short period planets (iHoward et al.ll2011 



Youdinll2011bl l. Characterization of the Kepler stars will thus powerfully constrain planetesimal 



formation models. 



5. PLANETESIMALS TO PLANETS 

Once planetesimals become larger than a few kilometers — potentially they are born much 
larger as discussed above — gravitationally focused collisions dominate growth into protoplanets. 
The size when "planetesimals" become "protoplanets" is vague. Although we use the terms inter- 
changeably, ~1000 km is a useful threshold. Depending on location and gas temperature, ~1000 km 
protoplanets are the smallest planets capable of binding disk gas into an atmosphere. 

We describe the accretion of solid protoplanets in §5.11 We start by discussing the processes 
that operate in a gas- free disk, including gravitationally focused collisions ( ^5.1.2p . velocity excita- 
tion ( ^5.1.3p and collisional fragmentation ( ^5.1.4p . Wh then describe planetesimal interactions with 
the gaseous disk ( §5.1.5p . Section [5.1.61 describes simulations of terrestrial planet formation that 
put these ingredients together. The accretion of a gaseous atmosphere ( ^5.2p affects planetesimal 
accretion ( §5.2.2p and transforms a planetary core into a gas giant ( §5.2.31 §5.2.4|) . We discuss nu- 
merical simulations combining the growth of giant planet cores and atmospheres in ^5.2.51 Finally, 
a young, massive gas disk might fragment directly into a gas giant or a brown dwarf. Section 15.31 
describes this formation channel and whether it might explain some exoplanets, especially massive 
giants at large radial distances. 



5.1. Growth of Solid Protoplanets 

Unlike planetesimal formation, it is easy to understand why planetesimals grow into larger 
protoplanets, even if the details are complicated. For the largest planetesimal in any region of the 
disk, collisions essentially always result in growth. Planetesimal velocities cannot be locally excited 
above the escape speed of the largest protoplanet. Consequently, the kinetic energy of collisions 
does not exceed the gravitational potential at the surface of the largest protoplanet. Collisions 
dissipate a fraction, sometime quite large, of the impact kinetic energy. Even if the impacting 
planetesimal shatters, growth is assured. 
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Collisions among smaller planetesimals, however, often lead to erosion or catastrophic frag- 
mentation. When an external, massive perturber stirs a belt of planetesimals, planetesimals collide 
at velocities larger than their escape velocity. These high velocity collisions tend to erode or com- 
pletely shatter planetesimals. The dust and changes in the planetesimal size distribution that result 
from these collisions are relevant for debris disks and for asteroids and Kuiper Belt objects. 

Unlike the planetesimal formation phase, aerodynamic drag no longer plays a starring role in 
protoplanet growth. However drag can still help regulate planetesimal velocities. The accretion of 
atmospheres (see ^5.2p also affects planetesimal capture. 

Deriving the precise evolution of a swarm of planetesimals is a complex numerical problem being 
attacked from several angles ( ^5.1.6p . However, we can develop a reasonabl y accurate picture of th e 
evolution with the "two groups approximation" reviewed in greater detail by Goldreich et al. ( 20041 ) . 
This approximation considers interactions between small, low mass planetesimals with mass and 
larger, more massive planetesimals with mass mi. The planetesimal masses ; = 47rr^j/9,/3 are 
related to their radii Vs^i and internal mass density, which we fix at p, = 2gcm~^ unless stated 
otherwise. Neighboring planetesimals have similar semimajor axes, a, and orbital frequencies, 0. 
Though detailed treatments need not make this approximation, we equate orbital eccentricities 
Cg^i and inclinations (in radians), but allow and ei to differ. The random velocities relative to 
a circular orbit are thus Vg^i ~ Cs^i^a and the vertical scale heights of the planetesimal disks are 
Hs,i ~ v^^ijfl. When the nature of the planetesimal is unspecified, we drop the s and I subscripts. 



5.1.1. Basic Length and Velocity Scales 
A useful scale for studying interactions of planetesimals and protoplanets is the Hill ( 18781 ) 



radius 

= {mT' " ■ <"> 

Planetesimals separated by < Rh are within the Hill sphere where their mutual gravitational 
attraction dominates the tidal gravity from the central star. While a mutual Hill radius can be 
defined, in practice it suffices to consider the more massive planetesimal. 
The size of a planetesimal in Hill units defines the parameter 

where is the mean mass density of the central star. Since 3p^/p, ~ 1, the parameter ip ~ R-k/a 
is roughly the angular size of the central star as observed from the planetesimal. The smallness of 
ip represents the fact that physical collisions are rare compared to gravitational scattering. 
At the Hill radius, the orbital speed about the protoplanet is the Hill velocity 

When planetesimal random velocities exceed vh, two-body encounters are "dispersion-dominated," 
negligibly affected by orbital shear. Random speeds below vh cause "shear-dominated" encounters 
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that involve (restricted) three-body dynamics. We describe below how the Hill velocity divides 
different accretion regimes. 

The outcome of a shear dominated encount er between planetesim als depends on the difference 
m semimajor axes 5R, relative to the Hill radius (jPetit HenoJlQSd ). When 6R < 1-2 Rh, the two 
planetesimals are deflected on a horseshoe orbit. More distant encounters with 6R > 2\/3Rh result 
in small angle scattering. For intermediate separations, 1-2 Rh ^ < 2\/^Rh, planetesimals enter 
the Hill sphere, experience chaotic deflections, and (if no collision occurs) leave the Hill sphere with 
typical relative velocity vh- 

In Hill units, the escape speed from the surface of a protoplanet, Vesc = [2Gm/r]^/^, is Vesc ~ 
vn/ip^^"^- Planetesimal velocities can be gravitationally excited up to the escape speed of the 
largest protoplanet. To estimate when a massive protoplanet might eject nearby planetesimals (or 
protoplanets) , we compare the escape velocity of the protoplanet to the orbital escape velocity, 
Vesc,i. ~ V2vK- When 

Vesc / m / a \m / 

^r'-'^iu;) (au) (d (^-^> 

exceeds unity, a planet of mass m can eject other nearby protoplanets. Terrestrial planets like Earth 
are too low mass and too close to the Sun to eject objects. The four Solar System giants can all 
eject pla netesimals; Jupiter is the most efficient at ejecting comets (and spacecraft) from the Solar 
System ( Fernandez &: Ip 19841 ). Aside from collisional grinding, eq. (j5-4p implies that planetesimal 
formation is more efficient closer to a star. 

The concepts of the Hill sphere and the Roche lobe are identical, though often used in different 
contexts. Both describe the region where the gravity of an object exceeds the tidal perturbation 
from its companion. Formally, both volumes are defined by the critical equipotential containing 
the Li and L2 Lagrange points. The Roche lobe is more distorted than a sphere when it describes 
binary stars that are similar in mass. The Roche radius or Roche limit describes the distance 
from a primary object at which the secondary becomes tidally disrupted and might form planetary 
rings. Aside from order unity corrections due to fluid effects or internal strength, the concept of 
individual vs. tidal gravity is again identical. To summarize, a secondary is at the Roche limit 
from the primary when it fills its own Hill sphere (or Roche Lobe). For an ensemble of very small 
planetesimals trying to become a much larger planetesimal, the Roche limit sets the critical density 
for gravitational collapse (eq. |4-13j l. 



5.1.2. Gravitationally Focused Collisions 



In this section, we describe the growth rates of large protoplanets (subscripted by I) accreting 
either large protoplanets or smaller planetesimals (unsubscripted). Gravitational focusing is the 
most important aspect of growth. Smaller random velocities for accreted planetesimals yield larger 
gravitational focusing factors and shorter growth times. We defer to later sectio ns the self consis- 
tent c al culation of planetesimal vel o cities and assume th e stand ard case, vi < Vg. iGreenberg et al 



(|l984l ). IWetherill fc StewartJ km± . iKenyon &: Bromleyl (j2008l ). and references in each paper de- 
scribe more detailed expressions for growth rates. 

We begin with the dispersion-dominated regime where v > vh,i- The mass accretion rate 
results from the usual isotropic expression as rhi = mnav. Adopting the surface mass density of 
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planetesimals, E, the number density of planetesimals, n, is 



mnv 



The cross section 



cr = 7r(n + rffcdisp > 



(5-5) 



(5-6) 



is the product of the geometric area and the gravitational focusing factor fcdisp- If the velocity of 
the incoming planetesimal at infinity is f > Vesc,h there is no gravitational focusing and fcdisp = 1- 
When V <C Vesc,h the speed on impact is roughly Vesc,i- Angular momentum conservation during a 
two body encounter sets the impact parameter for a grazing collision as riVesc,i/v- This expression 
yields fc. .disp ~ {'Vesc,l/v)'^ • Including energy conservation gives both cases simultaneously as 



fG,disp = ^ + PiVesc,l/vf 



(5-7) 



where /3 = 1 for a pur e two body interaction and ^2.7 accounts for anisotropic effects introduce d 
by orbital dynamics ( Greenzweig &: Lissauer 1990l : Spaute et al. 1991 : Wetherill fc Stewart 19931 ). 

Putting these results together and ignoring order unity coefficients gives the dispersion-dominated 
growth timescale mi /mi as 

tdisp ~ ^J!l ■ (5-8) 

^>'ijG,disp 

This result is just the geometric collision time over the focusing factor. 

For shear-dominated encou nters with v < y^-/, collision rates are affecte d by chaotic tra- 
jectories inside the Hill sphere ( Greenberg et al. 1991 : Pones &: Tremaine 19931 ). In this regime, 
planetesimal disks are thinner than the Hill radius, H ^ v/^l < Rh,i- Thus planetesimals enter 
the Hill sphere at the 2D mass accretion rate, rfiH ~ UR^if?. The probability, P, of a collision 
within the Hill sphere has two cases. Both use the maximum impact parameter for gravitationally 
focused collisions 6max ~ fiVesc,i/vH,i ~ 4'^^'^Rh,i- If the scale height of the disk is (relatively) thick 
with H V / Q > ftmaxi the collision probability P ~ 1)^^^/ {Rh^iH) is the ratio of the collision 
cross-section to the area of the accreting disk of planetesimals. For a thinner planetesimal disk, 
P ~ 

^max/ Rh,i is the ratio of the impact parameter to the Hill radius. 

Combining the mass flow rate through the Hill sphere with both limits of the collision proba- 
bility yields the shear dominated growth timescale, mi/{P ■ mn), as 



tshear 



shear 



(5-9) 



This timescale is again expressed as the product of the geometric collision time and a gravitational 
focusing factor 



fa. 



shear 



^ + ^3/2 



(5-10) 



Thus for V < iI^^^'^vh,i, gravitational focusing reaches its maximum value of fc ~ V"^''^) resulting 
in the fastest possible growth rate. Because inclination excitation is weak in the shear dominated 
regime, this fastest th in-disk accretion rate likely applies for all large protoplanets with vi < vh,i 



(jGoldreich et al.ll2004l ). Aside from this issue of anistropic velocities, the dispersion and shear 
dominated focusing factors match at f ~ vh,i with fc ~ 1/V'- 
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For numerical estimates of growth timescales we consider three cases. For the slow case we take 
V > Vesc,i and no gravitational focusing. For the intermediate case, we identify v ^ vh,i ~ '4^^^'^Vesc,i 
as the transition between shear- and dispersion-dominated. The fast case considers the maximum 
focusing factor fa ^ appropriate for v < 'ip^^'^VH,! ~ i^v^sci (and possibly for higher speeds 

when planetesimal i <^ e). Using eq. ()2-2p for the surface density of planetesimals, the growth times 
become 



^slow 
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(5-lla) 
(5-llb) 
(5-llc) 



These mass doubling times increase with protoplanet mass. Thus if gravitational focusing stays fixed 
or decreases (far from a certainty) these estimates also give the total accumulation time. These 
expressions omit the dependence on stellar mass and planet density for clarity but the growth 
times scale cx p^^^/M^^/^, pi^^ /M^,^/^ , and pi^^M^,^ for the three cases, respectively. Though 
stellar mass is not a major dynamical effect, it could correlate with disk mass or metallicity (here 
meaning planetesimal to gas ratio), normalized above by F and Zj-ei, respectively. Higher density 
protoplanets have smaller cross-sections and grow more slowly, but this effect becomes much less 
significant as gravitational focusing increases. 

Gravitational focusing dramatically speeds up the growth of protoplanets, especially in the 
outer disk. Without focusing, planets accumulate in tgiow ~ tens of Myr inside a few AU and more 
than 1 Gyr outside 5 AU. While a long growth time for terrestrial planets is acceptable, gas giants 
must form within a few Myr. Thus, formation of giant planet cores in the outer disk requires strong 
focusing, when the growth time for a 10 core at 50 AU, tfast ~ 3 Myr, is close to the lifetime 
of most gas disks. If strong f ocusing occur s, the formation of distant gas giants depends on the 
ability of cores to accrete gas (jRafikovl 120111 ) . 

Protoplanet accretion also depends on how the velocity distribution evolves. In dispersion- 
dominated gravitational focusing, the growth time tdisp oc l/rf, larger planetesimals grow faster 
than smaller ones. Although it is not the fastest regime, this "runaway" growth requires that 
gravitational focusing factors remain in the dispersion- dominated regime. With. Vp,^r/ vN ~ "0"^/^ oc 
a^/^, runaway growth persists longer in the outer disk ( Greenzweig Sz Lissauer 1990l ). 

When the largest protopanets enter the shear-dominated regime, runaway growth ends. For 
the thick disk case, t shear is independent of size. With either no focusing, t^iow ri, or the 
fastest (thin disk) shear-dominated accretion tf^st oc ri, smaller protoplanets grow faster and can 
catch up to the larger ones. In this "oligarchic" growth, many oligarchs compete to accrete small 
planetesimals, leading to an ensemble of oligarchs throughout the disk. 

As the largest protoplanets grow, they try to accrete all solid material in th eir vicinity. Tw o 
planets on circular orbits separated by a little more than 2\/3 Rh are stable (Gladman 19931 ). 
However , a fairly stable system with more planets requires larger separations, ~ BRh with B 
= 7-10 ( Lissauer 1987 : Kokubo &: Ida 19981 ). Planets that accrete all material within BRh are 
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"isolated." Setting rriiso = 27r'EaBRH leads to the isolation mass, 

With F = 1 and Zj-ei = 0.33, isolated objects in the terrestrial zone have masses comparable to 
Mercury and Mars. The MMSN has room for 30-50 isolated objects between the orbits of Mercury 
and Mars. Because their escape velocities are much smaller than their orbital velocities (eq. |5-4j ). 
isolated protoplanets eventually collide and merge to form Earth-mass planets (Fig. [SD. 

Outside the snow line, Z^ei = 0.78 at 5 AU yields an isolation mass of roughly 1 M^. As we 
show later, this mass is too small to bind the gas required for a gas giant. Increasing the mass of 
the MMSN {F ~ 5) increases the isolation mass to the 'typical' core mass of 10 M0 needed for 
a massive atmosphere. Thus, the MMSN is fine for the terrestrial planets, but it is not massive 
enough to allow formation of gas giants Simialr to Jupiter and Saturn. The extra mass required is 
consistent with observations of disks around the youngest stars (§2). 



5.1.3. Planetesimal Velocity Evolution 



As the previous section makes clear, the evolution of planetesimal velocities establishes the 
rate protoplanets accrete smaller planetesimals. Gravitational scattering is more common than 
physical collisions; thus, planetesimal velocities rapidly adjust as large protoplanets grow. 

Several processes modify the random velocities of planetesimals. The source of random kinetic 
energy is known as viscous stirring. This process uses planetesimal encounters — predominantly 
gravitational scattering — to extract energy from orbital shear. Dynamical friction redistributes 
kinetic energy among planetesimals of different masses, pushing them towards equipartition. Thus, 
small er (larger) planetesinials damp (excite) the random velocities of the la rger (smaller) planetesi- 
mals (jWetherih fc StewartJ[l989l : iKokubo fc Idall 199.4 IXenyon Sz Luull 19981 ). Ignoring ejections and 
gas drag, physical collisions are the only source of kinetic energy damping. Collis ional damping 
is especially effective for small planetesimals, r < 1-100 m, that collide frequently ( Ohtsukilll992l : 
Kenvon Luu 19981 ). Wh en collisions pro d uce s mall fragments that collide even more frequently, 
damping is very efficient. iGoldreich et al.l (120041 ) discuss order-of-magnitude derivations of these 
processes. As with accretion, behaviors vary between the dispersion- and shear-dominated regimes. 
It is common to refer to the excitation and damping of planetesimal velocities as "heating" and 
"cooling," respectively. 

The main goal of this introduction to velocity evolution is to show that planetesimals cannot 
be heated above — and can sometimes be cooled significantly below — the escape velocity of the 
large protoplanets. We focus on dispersion-dominated encounters to explain this result, which is 
crucial for ensuring the gravitationally focused collisions required to make planets on reasonable 
timescales. 

We first consider the simple case where all planetesimals have the same size. When v < Vesc, 
viscous stirring is dominated by gravitational scattering and occurs on the scattering timescale. 
This heating timesca l e is w ell approximated by the two body relaxation time from stellar dynamics 
( Binney Tremaine 20081 ). For the nav estimate of the gravitational scattering rate, we use 
eq. (j5-5p and compute the cross section (J ^ ^scatt froii^ tliG impact parametGr for strong gravita-tional 
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scattering, ftgcatt ~ Grajv^ . Together, these give the viscous stirring timescale ( Ida &: Making 19931 ) 



— ) , (5-13) 



where the constant C\ w 1/40 arises from a more detailed analysis ( Ohtsuki et al. 20021 ). and 



is similar to the Coulomb logarithm, In A, in stellar dynamics (and plasma physics). The final 
approximate expression in eq. ()5-13p facilitates comparison with the collision rates. 

The cooling rate is the gravitationally focused collision rate, which follows from eq. (I5-8P as 



Jesc 



Balancing the stirring and cooling rates implies v ^ v^sc- While the correct answer, the reasoning 
is incomplete. Gravitational focusing is weak for v > fg^c- In this regime, stirring and coUisional 
cooling rates are comparable. A slight imbalance in favor of heating could lead to runaway growth 
of V and an eventual collisional cascade. This runaway requires n early elastic physical co l lisions , 
as in the collision of two basketballs. In the idealized model of Goldreich &: Tremaine ( 19781 ). 



collisions among planetesimals with coefficients of restitution > 0.63 (comparable to a baseball, but 
smaller than a basketball or a table tennis ball) bounce often enough to lead t o runaway heating 



Coefficients of restitution for planetesimals are probably much smaller than 0.5 (jPorco et al.ll2008l ): 
the velocity runaway is unlikely. Similarly sized planetesimals will excite random velocities to the 
surface escape speed, v ~ Vesc- 

Returning to the two groups approximation, we consider stirring of small planetesimals by 
large protoplanets. Dynamical friction ensures Vs > vi (confirmed below); planetesimals dominate 
the encounter speed. The stirring of small planetesimals by larger ones then occurs on a timescale 

, P,ri ( Vs \^ ,^ 

t-stir.disp ~ ■c:^ ^ I I ■ 1^0- iOj 

luiil \Vesc,lJ 

Comparison with eq. (|5-13p shows that large planetesimals dominate the stirring of small planetes- 
imals when XJiTJii > Us'^s- Initially, Us > Ui; small planetesimals contain enough mass to affect 
growth. To dominate stirring, however, large planetesimals can contain a minority of the mass. 

Due to stronger stirring by large protoplanets, small planetesimals are excited to Vg > Vesc,s- 
At these speeds, collisions between small planetesimals generally cause collisional fragmentation or 
erosion. The resulting smaller planetesimals then collisionally cool more efficiently. Even without 
this extra cooling, gravitational focusing arises. When Vg > Vesc,s: small planetesimals cool by 
colliding with other small planetesimals on the geometric timescale, tcooi ~ P»fs/{Us^)- Balancing 
these heating and cooling rates gives 

Vs ^ "F -] '"esc,l ■ (5-16) 

\^s n ) 

When i7s > ^ ^esc,/; small planetesimal accretion is strongly gravitationally focused. 

We now consider whether the growth of large protoplanets is dominated by the accretion of 
small planetesimals or other large protoplanets. Planetesimals with the larger product of surface 
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density and gravitational focusing, Ufa, drive the fastest growth (eq. |5-8H5-9| ). A balance of 
viscous self-stirring and cooling by dynamical friction against small planetesimals then sets the 
velocity dispers ion of large protoplan ets. For dispersion-dominated encounters this balance gives 



(for details, see iGoldreich et al.ll2004l ) 



(5-17) 

Since Sfcdisp oc S /v"^ , small planetesimals contribute more to the growth of large protoplanets, 
by a factor [Es/I^iY/'^ > 1. 

This introduction only begins to touch on the complexities of planetesimal velocity evolution. 
However even these simple considerations show that gravitationally focused accretion of small 
planetesimals by large protoplanets is likely. Earlier, we explained that collisional erosion plays a 
key role in cooling small planetesimals to Vg > Vesc,s- Now, we turn to even more violent encounters, 
with Vs S> Vescs, which lead to catastrophic disruption. 



5.1.4. Fragmentation 



As large planetesimals grow, they stir up the velocities of smaller planetesimals to the disrup- 
tion velocity. Instead of mergers, collisions then yield smaller planetesimals and debris. Continued 
disru ptive collisions lead to a collisional cascade, where leftover planetesimals are slowly ground to 
dust (jPohnanvil 1 19691 : 1 Williams k Wethe"riilll994l : Ixobavashi fc Tanakalboiol '). Radiation pressure 
from the central star ejects dust grains with r < 1-10 //m; Poyii t ing-Robertson drag pulls large r 
grains into the central star ( Burns et al. 19791 : Artymowicz 1988 : Takeuchi &: Artymowiczl I2OO1I ). 
Eventually, small planetesimals are accreted by the large planetesimals or ground to dust. 

To understand the origin of the collisional cascade, we consider the outcome of a head-on 
collision between two identical planetesimals. During the impact, some kinetic energy heats up 
the planetesimals; the rest goes into the internal energy of material in the planetesimals. When 
the impact energy is small, the extra internal energy is small compared to the binding energy of 
either planetesimal; the two objects merge into a single, larger planetesimal. When the impact 
energy is larger than the binding energy, the collision shatters the planetesimals into a few smaller 
planetesimals and a lot of dust. 



Estimating the binding energy of pl anetesima ls relies on two approa c hes (IDavis et al 



1985 



Housen fc Holsapplel [l99ol Tl999l:lHolsaPDlc 1994 iBenz Asphaud Il999l : iLeinhardt etail l200^ 
Leinhardt fc Stewart 20091 ). Sophisticated collision experiments yield the internal strengths of 
small rocky and icy objects, r < 10-100 cm. Theoretical investigations derive the strength from 
analytic or numerical models of the crystalline structure and the equation-of-state of the material. 
In both cases, investigators derive the energy Q*^ required to remove half of the combined mass of 
two colliding planetesimals and eject this mass to infinity. Although more sophisticated approaches 
include the impact velocity in Q^, we focus on a simpler expression that depends only on radius. 



J9 



Here Qhfs^ is the bulk (tensile) component of the binding energy and pQgt 
ponent of the binding energy. 



(5-18) 

is the gravity com- 
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Laboratory experi ments and detailed numerical collision simula,tions yield a broad range of 
results for Q*p (Fig. [2 iHousen k Holsappld Il990l : iBenz Asphaud Il999l : iHolsapple et all l2002l : 



Leinhardt et al. 20081 ) . In the strength regime at small sizes, the binding energy of a planetesimal 
depends on the number of flaws - cracks, fissures, etc - in the material. Larger planetesimals have 
more flaws and smaller strengths. In the gravity regime at large sizes, the binding energy depends 
on the internal pressure. Larger planetesimals have larger internal pressures and larger strengths. 
The lower density and weaker crystalline structure of ice leads to smaller strengths than basalts 
and other rocks. 

Mod els for the breakup of c omet Shoemaker-Levy 9 suggest a smaller component of the bulk 
strength ( Asphaug &: Benz 19961 ). implying small disruption energies for small planetesimals (Fig. 
[71 "Rubble Pile"). A low strength is consistent with numerical simulations of "rubble piles", 
structures with countless flaws held loosely together. This structure probably results after icy or 
rocky planetesimals suffer numerous impacts which disrupt the internal structure (removing most 
of the tensile component of the binding energy) but do not destroy the object. 

The collisional cascade begins when the impact energy of colliding small planetesimals equals 
Q*ij. Because the random velocities of small planetesimals equal the escape velocities of large 
planetesimals, the impact energy depends only on the mass of a large planetesimal. Equating 
this energy to Qd* allows us to derive the "disruption mass," the mass of a large planetesimal at 
the onset of the collisional cascade. With Vesc,s ^ ^^esc,Z! the impact energy per unit mass in the 
center-of-mass frame is roughly 
disruption mass: 



Setting this energy equal to Q*^, we solve for the 



rrid 



1/2 



8Q 



D 



G 



3/2 



3.5 X lO^V. 



6^-1/2 



107 



erg g- 



3/2 



Ma 



(5-19) 



When small planetesimals have sizes exceeding ~ 1 km, is fairly independent of their compo- 
sition. For typical ~ lO'^ ~ 10^ erg g~^, the disruption mass is roughly 0.003-3 Pluto masses. 
Collisional cascades begin well before planets reach their final masses. 

Once disrup tion commences, the final mass o f a planet depends on the times cale for the col- 
lisional cascade ( Kenyon fc Bromley 2004a . 20081 : Leinhardt Sz Richardson 20051 ). If disruptive 
collisions produce dust grains much faster than planets accrete planetesimals, planets cannot grow 
much larger than the disruption radius and have a maximum mass mi^rnax ~ w-d- However, if 
planets accrete grains and leftover planetesimals effectively, planets reach the isola tion mass before 
collis ions and radiation pressure remove material from the disk {mi^max ~ nT-iso] iGoldreich et al 



20041 ). 

In a gas- free environment, larger in the inner disk enables planets to accrete much of the 
debris before destructive collisions and radiative processes remove it. Rocky planet masses then 
approach the isolation mass. In the outer disk, smaller planets cannot accrete debris before it is 
lost. Icy planets cannot grow much larger than rrid- 



5.1.5. Planetesimal Accretion with Gas Damping 

Gas slows the random velocities of smaller planetesimals. Larger protoplanets are less affected 
by drag and are damped by dynamical friction. The drag force, Fd, of eq. ()4-4p damps the kinetic 
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energy of planetesimals (now with size r^, not the size s of dust, pebbles, and boulders) at a rate 
tgas = Vs{dvs/dt)~^ = Vs{F£)/ms)~^. With pg ~ Ugf2/cs, where Cg is the gas sound speed, 



1 p,r= Cs 



To understand the impact of damping, we consider an ensemble of small planetesimals stirring 
themselves (eq. |5-13| ). Without gas, small planetesimals excite their velocities to Vs ~ Vesc,s- If 
tgas < tcoii ~ PrnTs/i^s^), drag exceeds collisions as the dominant cooling mechanism. This switch 
happens when 

To be damped by gas drag, planetesimals must exceed this minimum size. This somewhat counter- 
intuitive result (drag is often more significant for smaller bodies) arises from (i) non-linear drag laws 
and (ii) a velocity scale, Vesc,s, that increases with size. The numerical value of the size threshold 
decreases if Z^^i is reduced due to an inefficiency of turning dust into planetesimals. 

When eq. (j5-2ip holds, the stronger damping of self-stirred planetesimals ensures Vg < Vesc,s- 
Collisions are gravitationally focused and runaway growth begins earlier than in a disk without gas. 

As growth proceeds, larger protoplanets dominate the stirring of smaller planetesimals. With 
stronger stirring, smaller and smaller planetesimals are damped by non-linear gas drag instead of 
collisions. Even if collisions are initially more significant, gas drag becomes the dominant coolant 
as growth proceeds. To compute the random speeds of small planetesimals, we assume dispersion- 
dominated encounters and balance the heating of eq. (I5-15P with the cooling of eq. ()5-20p to get 

1 r.U,V/' , , 

. (5-22) 



\Cd n 

With more mass in small planetesimals < <C X'g, gravitational focusing, fcdisp ~ {vesc,i/vs)'^, 
becomes strong. Accretion of small planetesimals can become shear-dominated. As described in 
^5.1.2| runaway growth transitions to oligarchy (but does not slow down) in the transition to shear- 
dominated accretion. A self-consistent analysis of these processes is facilitated by the numerical 
calculations summarized in ^5.1.61 

Details aside, the large oligarchs stir small planetesimals past their escape speed and up to the 
disruption velocity (§5.1.3). Disruptive collisions among small planetesimals produce a collisional 
cascade, which grinds planetesimals into smaller and smaller objects. Without gas, planetesimals 
are ground into small dust grains which are dragged into the star by Poynting-Robertson drag 
or ejected from the planetary system by radiation pressure. With gas damping, the collisional 
cascade halts at some intermediate size (~ 10 cm to 10 m), depending on factors such as the mass 
of the oligarchs, gas density, material strength, and orbital distance. The damped velocities are 
then slow enough that the oligarchs accrete these small rocks rapidly. This rapid ac cretion enables 



oligarchs to reach the is olation mass on short timescales, even in the outer disk (jRafikovl 120041 : 
Kenyon &: Bromley 20091 ) . 



Collisional grinding a set of small planetesimals into small dust grains requires a very depleted 
gas disk. For the cascade to proceed down to 1-10 ^m particles, these grains must decouple from the 
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gas (on an orbital timescale, we assume for simplicity). Epstein drag applies for low gas densities 
(and long gas mean free paths). From eq. (j4-3a|) . particles with sizes s > ^g/ p» decouple from the 
gas. From eq. (|2-ip the depletion factor (relative to the MMSN) required to avoid entrainment is 

This low mass disk may not last very long. Nevertheless, current observational limits only constrain 
gas surface densities in debris disks to F < 1% of the MMSN, not yet sufficient to asses the 
dynamical significance of gas. ALMA should should place much tighter constraints (see chapter by 
Moro-Martin). 



5.1.6. Numerical Simulations of Low Mass Planet Formation 



Analytic estimates provide a good understanding of each piece of the planet formation process. 
However, putting the whole set of processes into a coherent theory requires numerical calculations. 
Clusters of computers can now finish an end-to-end calculation in a reasonable amount of time. 
Several groups are building towards this simulation, but no complete calculation exists. 

Constructing numerical simulations of planet formation involves identifying and solving a set 
of coupled differential equations which describe the evolution of the gaseous disk and the masses 
and orbital properties of solid objects. Selecting the proper approach depends on the nature of 
the problem. Hydrodynamics codes address t he evolution of the gaseous disk and how planet s 



accrete material and migrate within the disk ([D'Angelo et al.l l2003l : iNelson fc Papaloizoul |2004| ) 



Smooth particle hydrodynamics allows detailed solutions to outco mes of binary collis ions between 
large protoplanets (e.g., Earth-Moon and Pluto-Charon formation: ICanupll2008l . I2OIII ). Solving the 



coagulation equation with a fragmentation algorithm yi elds the inass and time evolution of solid 
particles rangiiig in s i ze from 1 ^m up to rq ug hly 1000 km (ISafronovlE969l : IWetherill StewartlEooa : 
Kenvon fc Luulll999l : iBirnstiel et ai1l2010l l. To treat the d ynamical evolution of large planets, A^- 
body treatments provide accurate and often fast solutions (IChainbers &: Wetheril]|[l998l : IChamberi 
2OO1I : iRavmond et"aDl2004l : iNaeasawa eraDl2005l : iKokubo et ahlbood l. 

Most investigations of terrestrial planet formation employ a coagulation code or an A^-body 
code. An A^-body code cannot possibly follow the trajectories of the > 10^^ small planetesimals 
expected in a MMSN. Coagulation models, which treat planetesimals as a statistical ensemble of 
objects with a distribution of e and z, can s olve for the time evolution of their masses and orbits 
throughout runaway and oligarchic growth ( Wetherill &: Stewart 19931 ). Once most of the solid 
mass is in a few protoplanets, the statistical approach fails. A^-body codes can then follow the 
evolution during the late stages of oligarchic growth and throughout chaotic growth. 

S everal hybrid codes combi i ie aspects of both approaches (ISoaute et al I1991I1 [Weidenschilling et al 



19971 : iBromlev KenvonI l2006l : ICharnoz fc Morbidellil 120071: iRavmond et al 



evolution of the gaseous disk together with the solids, iBromley &: Kenyon 



2011 



(2011 



). To follow the 
solve the radial 

diffusion equation for the gaseous disk (eq. |3-4j ) and employ a merged coagulation -|- A^-body code 
for the solids. In these treatments, the coagulation code follows solids with masses smaller than the 
promotion mass, rupro', the A^-body code tracks protoplanets with m > rupro- Co mparisons with 



other simulations and with analytic theor y provide tests of these techniques (e.g., iKenvon &: Luu 



19981 : lFraseill2009l : iMorbidelh etaDl2009bl ) 
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To illustrate the formation process, we summarize results for several calculations of terrestrial 
planets and gas giant cores. Because this aspect of this field is growing so rapidly, we focus on a 
few simple examples. 



^Coagulation codes begin with an ensemble of planetesimals in place at t = ([Kenvon &: Bromley 

Planetesimals are placed in concentric annuli according to a fixed initial surface density re- 



lation. These planetesimals often have a single size of 1-100 km; sometimes calculations begin with 
a distribution of sizes. Because dynamical friction efficiently damps the velocities of the largest 
planetesimals, planets grow faster in calculations with a size distribution of planetesimals. Starting 
with an ensemble of small planetesimals leads to faster growth than an ensemble of large planetes- 
imals. The initial surface density sets the growth time. Planets grow faster in more massive disks. 
In many calculations, the planetesimals evolve in a gaseous disk which also evolves in time; the 
disk evolution may be proscribed in advanced or calculated along with the planetesimals. 

Fig. [8] shows the evolution of oligarchs in an evolutionary sequence starting with an ensemble 
of 1 km planetesimals at 0.4-2 AU. Following a short runaway growth phase, protoplanets with 
m > rripro appear in a wave that propagates out through the planetesimal grid. As these oligarchs 
continue to accrete planetesimals, dynamical friction maintains their circular orbits and they evolve 
into "isolated" protoplanets. Eventually, large oligarchs start to interact dynamically at the inner 
edge of the grid; a wave of chaotic interactions then moves out through the disk until all oligarchs 
interact dynamically. Once a few large oligarchs contain most of the mass in the system, dynamical 
friction between the oligarchs and a few leftover planetesimals starts to circularize their orbits. 
This process excites the lower mass oligarchs and leftover planetesimals, which are slowly accreted 
by the largest oligarchs. At the end of the calculation, the masses, semimajor axes, and orbital 
eccentricities of stable planets are similar to those of the terrestrial planets in the Solar System. 

Comparisons between the results of hybrid and A^-body calculations show the importance of 
including planetesimals in the evolution. Both approaches produce a few terrestrial mass planets in 
roughly circular orbits. Because dynamical friction between leftover planetesimals and the largest 
oligarchs is significant, hybrid calculations produce planets with more circular orbits than tradi- 
tional A^-body calculations. In most hybrid calculations, lower mass planets have more eccentric 
orbits than the most massive planets, as observed in the Solar System. In both approaches, the 
final masses of the planets grow with the initial surface density; the number of planets is inversely 
proportional to the initial surface density of solid planetesimals. However, the overall evolution is 
faster in hybrid calculations: oligarchs start to interact earlier and produce massive planets faster. 

In hybrid calculations, the isolation mass and the number of oligarchs are more important as 
local quantities than as global quantities. As waves of runaway, oligarchic, and chaotic growth 
propagate from the inner disk to the outer disk, protoplanets growing in the inner disk become 
isolated at different times compared to protoplanets growing in the outer disk. Thus, the isolation 
mass in hybrid models is a function of heliocentric distance, initial surface density, and time^ which 
differs from the classical definition (eq. |5-12j ). 

During oligarchic growth of the simulation in Fig. [HI viscous stirring excites leftover planetes- 
imals to the disruption velocity. A series of separate simulations demonstrates that the collisional 
cascade produces copious amounts of dust, which absorb and scatter radiation from the central 
star. Following the growth of protoplanets, the cascade begins at the inner edge of the disk and 
moves outward. For calculations with a solar-type central star, it takes ~ 0.1 Myr for dust to form 
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throughout the terrestrial zone (0.4-2 AU). The timescale is ~ 1 Myr for the terrestrial zone of an 
A-type star (3-20 AU). As the collisional cascade proceeds, protoplanets impose structure on the 
disk (Fig. [9l left panel). Bright rings form along the orbits of growing protoplanets; dark bands 
indicate where a large protoplanet has swept up dust along its orbit. In some calculations, the dark 
bands are shad ows, where optically thick dust in the inner disk prevents starligh t from shining on 
the outer disk (Groean et al.ll200ll : iKenvon fc Bromlevlliooial : IPurda et al]l2004l ). 

In the terrestrial zones of A-type and G-type stars, the dust emits mostly at mid-IR wave- 
lengths. In calculations with G-type central stars, formation of a few lunar mass objects at 0.4-0.5 
AU leads to copious dust production in a few thousand years (Fig. [U right panel) . As protoplanets 
form farther out in the disk, the disk becomes optically thick and the mid-IR excess saturates. Once 
the orbits of oligarchs start to overlap (~ 1 Myr), the largest objects sweep the disk clear of small 
planetesimals. The mid-IR excess fades. During this decline, occasional larg e collisions generate 
large clouds of debris that produce remarkable spikes in the mid-IR excess (jKenyon &: Bromley 
20021 . 12OO5I I. 

In A-type stars, the terrestrial zone lies at greater distances than in G-type stars. Thus, debris 
formation in calculations with A-type stars begins later and lasts longer than in models with G-type 
stars (Fig. [9l right panel). Because the disks in A-type stars contain more mass, they produce 
larger mid-IR excesses. At later times, individual collisions play a smaller role, which leads to a 
smoother evolution in the mid-IR excess with time. Although the statistics for G-type stars is 
incomplete, current observations suggest that mid-IR excesses are larger for A-type stars than for 
G-type stars (see Chapter by Moro-Martin) . 

Collisional cascades and debris disk formation may impact the final masses of terrestrial plan- 
ets. Throughout oligarchic growth, roughly ~ 25% to 50% of the initial mass in planetesimals is 
converted into debris. For solar-type stars, the disk is optically thick, so oligarchs probably accrete 
the debris before some combination of gas drag, Poynting-Robertson drag, and radiation pressure 
remove it. In the disks of A-type stars, the debris is more optically thin. Thus, these systems may 
form lower mass planets per unit surface density than disks surrounding less massive stars. Both 
of these assertions require tests with detailed numerical calculations. 



5.2. Accretion of Atmospheres 

Protoplanetary atmospheres have a rather different character than the mature planetary atmo- 
spheres of Solar System planets and exoplanets. The crucial distinction is that protoplanets orbit 
within a gas disk. The disk supplies the atmosphere's gas and provides an external binding pres- 
sure until the planet opens a clean gap. The protoplanetary atmospheres inherits its composition 
from the disk. However, the fraction of heavy elements that wind up in the core - versus the dust 
and ablated planetesimals that remain in the atmosphere - is a key uncertainty. This uncertainty 
crucially affects the mean molecular weight, fi, and opacity, k, of the atmosphere. 



5.2.1. Static Protoplanet Atmospheres 

As protoplanets grow, they become massive enough to bind a gaseous atmosphere. The atmo- 
sphere is significantly denser than the surrounding disk gas when the core's gravitational binding 
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energy exceeds the thermal energy of the gas. For a solid core of mass nic and radius rc, this occurs 
when rB > where the Bondi radius 

rB = Gm/cl , (5-24) 

where total planet mass m = rric + ma, including the gravitationally bound atmosphere's mass, ma- 
Equivalently a core with an atmosphere must exceed the Bondi mass 

3 / ^ X -Q/14 



c: 

rric > niB = 



where we use the gas temperature in an irradiated disk (eq. |3-29| ) and normalize the stellar 
luminosity as 1^, = L^/L© and the gas mean molecular weight as /i = ^/(2.4 mn)- 

As the core mass increases beyond rriB the atmosphere becomes denser and more massive. 
The outer boundary of the atmosphere, rout = ^^^{fB, Rh), is set by the Bondi radius until the 
atmosphere fills the Hill sphere. This transition occurs for massive protoplanets, as > 
requires 



m > mt 



r.3 . r, . 6/7 ;3/8 



rans 



^3 — -=zTr. Me- (5-26) 



Comparing to the disk's gas scale- height Hg = c^/il, the criteria rs > Hg and Rh > Hg also repro- 
duce eq. (I5-26P within order unity constants (mtrans increases by 3'^/^ for the Rh > Hg criterion). 
Thus when m > m-trans the protoplanet is no longer uniformly embedded in the disk midplane. It 
can feel the top and bottom of the disk and start to open a gap (see chapter by Morbidelli). Outside 
~ 5 AU, the core accretion instability generally occurs for m < m-trans- Thus the 3D structure of 
the gas disk can usually be ignored when describing the onset of core accretion ( §5.2.3p but not its 
final evolution ( g5.2.4p . 

The structure of an proplanetary atmosphere obeys the equations (which also govern stellar 
structure) of hydrostatic balance 

dP Gm 
dr 

mass conservation 



-2-P = -P9 , (5-27) 



^ = Anr^p , (5-28) 
dr 



and energy transport by optically thick, r ~ nP/g > 1, radiative diffusion, 

16o-sbT^ dT _ L 
3kp dr Anr"^ 



(5-29) 



Wherever radiative diffusion would satisfy the Schwarzchild criterion, dlnT/dlnP > Vad, the 
energy transport becomes convective. The adiabatic index Vad = 2/7 for a n ideal diatomic gas , 
but in general must be d etermined from detailed equation of state calculations (jSaumon et al. Il99d; 



Saumon &: GuillotI |2004| ) . Because convective transport is efficient, the temperature profile follows 
an adiabat, T oc P^^d instead of eq. (j5-29j) in convective regions. For the ideal gas equation of 
state, P = pTZT, with the (specific) gas constant TZ = ks/ pmn- 

The masses of stable atmospheres (that do not undergo the core accretion instability) place 
interesting constraints on planet formation. In the terrestrial zone, typical isolation masses (eq. 
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|5-12j ) are roughly 0.1 M^. For almost any accretion time, these planets have stable atmospheres 
with masses much smaller than the planet's mass. Icy planets formed at tens of AU, however, have 
much larger isolation masses of several M®, and can support much more massive atmospheres. 
In the Solar System, the dichotomy between terrestrial planets with thin atmospheres and icy 
planets with massive atmospheres is consistent with our estimates. Once we have a large sample 
of rocky/icy exoplanets with well characterized atmospheres, we will see if the same dichotomy 
persists. 



5.2.2. Enhanced Planetesimal Accretion 
For low mass protoplanets with rout = the size of the atmosphere relative to the core is 

For Mars-mass planets and larger, the radius of the atmosphere is 10 or more times larger than the 
radius of the core. Extended atmospheres can significantly enhance planetesimal accretion. 

When a planetesimal encounters the atmosphere of a protoplanet, it experiences enhanced gas 
drag. Capture results if the planetesimal loses enough orbital energy. A thin atmosphere has little 
impact on very large planetesimals; the collisional cross-section is still nr'^fg. For sufficiently small 
planetesimals, any encounter with the atmosphere allows the planet to capture the planetesimal; 
the collisional cross-section is then Trr^^^/g. For intermediate siz es, the effective cross- section lies 
somewhere between vrr^/^ and T^r^uifg- To address this regime, llnaba fc Iko^ tooA ) define an 



"enhanced radius" r^, where the collisional cross-section is Tirlfg. In this approach, r^ = r for 
accreting very large planetesimals and r^ = rout for accreting very small planetesimals. However 
very small rocks and dust grains will be too tightly coupled to the gas to accrete. 

Compared to the estimates for dispersion- and shear- dominated growt h in §5.1, atmospheres 
can enhance accretion rates by 1-2 orders of magnitude ( Chambers 20081 ). When leftover small 



planetesimals have typical radii of 0.1-10 km, an isolated terrestrial planet with a thin atmosphere 
has a small radius enhancement, rg w 1 — 2. Thus, these isolated objects never experience rapid 
growth from an enhanced radius. Isolated icy planets are more massive and support massive 
atmospheres. These objects have r^ ~ 10 for accreting 0.1-10 km planetesimals and rg ~ 3 for 
accreting 100+ km planetesimals. Once they develop atmospheres, icy isolated objects rapidly 
sweep up any leftover planetesimals along their orbits. 



5.2.3. The Core Accretion Instability 

Low mass protoplanets (near m^) have low mass, optically thin atmospheres. More massive 
cores bind thicker atmospheres, which trap the heat of planetesimal accretion. As the heat is 
radiated away, the atmosphere becomes denser and more massive. When the atmosphere's mass 
exceed s roughly the core mass, th e atmosphere cannot maintain hydrostatic equilibrium and col- 
lapses (|Harrislll978l : lMizunolll98ol l. This collapse — referred to as the "core accretion instability" 
— leads to rapid gas accretion and the birth of a gas giant. The critical core mass (sometimes 
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called the "crossover mass" because collapse occurs when the core and atmosphere masses rou ghly 
match) for this instability is often, but not always, ~ 10 Me (|lkoma et alJbood : iRafikovlbood ). 



Near the crtical core mass, real protoplanetary atmospheres are convective in the interior and 
(for l ow enough pla netesimal accretion rates) radiative in the exterior region that match es onto the 
disk (|Rafikovll2006l ^. However an illustrative, and historically important, calculation by IStevensonI 



(|l982l ) demonstrates the essential features of the core accretion instability by (incorrectly) assuming 
the atmosphere is completely radiative with constant opacity n. We summarize this calculation 
before comparing it to more detailed computations. 

Following Stevenson ( 19821 ). we calculate the structure and mass in the atmosphere by keep- 



ing the mass, m, constant in the hydrostatic balance equation (|5-27p . i.e. neglecting the detailed 
variation from m = rric at the core to m = rric + rua at the top. Equations (I5-27P and (15-29P then 
give (temporarily omitting order unity coefficients for clarity) 

T^dT/dP ~ KL/{(TsBGm) . (5-31) 

To integrate this equation, we assume (correctly) that P and T are significantly higher at the base of 
the atmosphere than in the disk. We further keep L constant, appropriate if accreted planetesimals 
release their kinetic energy at the core's surface. This assumption yields L = GnicTn/rc, where rh is 
the planetesimal accretion rate. The slowing of planetesimals as they fall through the atmosphere 
is a minor correction included in detailed numerical models. 

With these approximations, eq. ()5-3ip integrates to a simple T — P profile 

T-(^^pY' . (5-32) 

To obtain the density profile, we use the fact that for a barotropic relation P cx r^°° (Voo = 1/4 
for our example of a constant opacity) , the hydrostatic balance equation for an ideal gas integrates 
to 

r = Voo^ (5-33) 
in the atmospheric interior. The density profile follows as 

fJsB / Gm\^ 1 



More generally, p (x i/Voo ^ -^rhich is left as an exercise. 

The atmosphere's mass follows by integrating eq. (j5-28p from rc to rout- 



(^SB f Gm\'^ ^ ^(JsBX 



4 



with X = lii(ro„j/rc) and order unity coefficients reinstated (despite the overall approximate nature 
of the calculation). The final expression relates the protoplanet luminosity to the (current) core 
growth timescale, tacc = nic/m = Grn^/ {rcL). Setting m = rric we can numerically evaluate 

-a-9.4(-^)'^'^^ Me, (5-36) 
V 10 Me / Ki Gyr 
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with X = 6.7 (from eq. |5-30| ) and ki = k/(1 cm^ g^^)- For the chosen parameters, we are near the 
crossover mass, rua ~ rric ~ 10 Mq. Fig. [TU] shows the behavior of this simple atmosphere model. 
We emphasize that the numerical values of this simple model are only meant to be illustrative. For 
instance the core must actually accrete in tacc ^ 10 Myr <C Gyr. 

The existence of an instability arises from the non-linearities in ma- Expressing 

m = rric + k — —p- (5-37) 

rUc' 

where k incorporates all the constants in eq. (j5-35p . we can show that beyond a critical core mass 
the total mass (unphysically) declines as nic increases. Using calculus, the turnoveiH occurs where 
drric/dm = at m = mc'^^(4/c)~^/^ and nic = 0.19A;~^/^. 

This simple derivation at least qualitatively explains many features of more detailed core 
accretion models. The strong dependence of the atmospheric mass on u is suppor ted by studies 
showing that envelope pollution lowers the critical core mass dHori Ikoinall201lh . The opacity 



is very sensitive to the amount and sizes of dust g rains ( Pollack et al. 19851 ). One popular way to 



speed up core accretion is to reduce the opacity (jHubickyj et al.ll2005l ). The "correct" choice of 
opacity likely varies between planets and is poorly constrained. It is unclear how much dust from 
ablated planetesimals will remain in the atmosphere. Small grains both contribute significantly to 
opacity and settle slowly. 

High planetesimal accretion rates increase the critical core mass. For very high accretion 



rates , especially in the inner disk, the protoplanet atmosphere will be fully convective (iRafikov 



20061 ). The atmosphere then matches onto the same adiabat (constant entropy curve) as the disk 



gas and thus has the lowest possible mass. In this case the formation of gas giants is quite unlikely. 
If planetesimal accretion stops (or becomes suitably small), the relevant luminosity comes from 
the Kelvin-Helmolz contraction of the atmosphere. Models that omit planetesimal accretion (but 
correctly compute contract ion) thus provide a meaningful lower limit on the critical core mass 



(jPapaloizou Nelsonll2005l ). 

It has long been postulated that the cores of Solar System giants might correspond to an 
isolation mass, ^ 5- 20 M^, at 5-10 AU, even though this would require a massive planetesimal disk 



(jPollack et al.lll99fil ). To understand why this assumption is reasonable, we turn to Fig. [TOl Initially 



a low mass planet accretes planetesimals rapidly. With a short accretion time, the puffy atmosphere 
is well below the crossover mass. As the core grows in mass, two effects bring the atmosphere closer 
to instability, (i) extra compression of the gas and (ii) fewer and fewer planetesimals to accrete 
and heat the atmosphere. As the atmosphere cools, the critical core mass drops until it reaches 
the actual core mass. This cooling is most likely to happen after most planetesimals have been 
accreted, i.e. at the isolation mass. 



Our values differ slightly from IStevensonl l| 19821 ) because we assume constant accretion time instead of constant 
mass accretion rate. Further the 3/4 exponent in his eq. (15) is a typo that should be 3/7. 
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5.2.4- Direct Accretion of Disk Gas (and How it Stops) 



When the envelope collapses, the planet starts to accrete gas directly from the protostellar 
disk. This dynamical process is not amenable to a stellar structure calculation. Direct accretion of 
gas is similar to accretion of planetesimals by oligarchs; the accretion rate is roughly the product 
of the planet's cross-section (or impact parameter in the 2D limit), the local gas density (or surface 
density), and the encounter velocity, which is dominated by Keplerian shear for circular orbits. 

Shortly after the core accretion instability, the planet crosses the transition mass (eq. |5-26j ) 
and the relevant accretion radius is the Hill radius, Rh. Since the transition mass also corresponds 
to the Hill radius exceedin g the gas s cale height (see discussion after eq. |5-26| ) accretion is not at 
the classic Bondi rate (e.g. IShulll992l ). Instead the two dimensional mass accretion rate applies, 



rup 



E^QRjj ~ 10^ 



F 



m 



1/6 V 60 Mf, 



2/3 



5 AU 



Me 
Myr 



(5-38) 



At 5 X 10 ^Mq yr i, this rate exceeds the accretion rate onto most pre- main sequence T Tauri 
stars! 

Clearly, something must stop this influx of gas. The gaseous isolation mass is one natural 
stopping point. Applying eq. (|5-12p to the gas disk gives 



mi. 



{2^^BE^^, ^ ^F^ 
(3M^)i/2 ^ 



a \3/4 
5 AU/ 



Mj 



(5-39) 



While this result appears plausible for Jupiter, most disk models for Jup iter's core require F > 5, 
requiring a second mechanism to halt gas accretion ( Lissauer et al.ll2009l ). 



Opening a gap in the disk — which begins when eq. (j5-26p is satisfied — can slow down 
accretion so that the disk dissipates before the planet reaches the gaseous isolation mass. Since 
disk lifetimes are at least at least 1-3 Myr, accretion times must be at least this slow to halt growth 
(without invoking a fine-tuning of core accretion and disk dissipation timescale s). Only a wide and 
relat ively clean gap can slow accretion enough to explain final planet masses (jD'Angelo &: Lubow 
20081 ). The effective viscosity of t he disk must be low f or a clean gap, which is an especially strong 
concern for self-gravitating disks (iKratter et al.ll2010bl . §5.3p . 



5.2.5. Numerical Simulations of Gas Giant Planet Formation 
To conclude this section. Fig. [TT] illustrates the evolution of the semimajo r axes of icy and gas 



giant planets from one simulation of material outside the terrestrial zone (e.g.. iBromlev &: Kenvon 



20 111 ). The calculation begins with a single 1000 km planetesimal and an ensemble of ~ 1 cm 
planetesimals in each of 96 annuli from 3 AU to 30 AU. Because the system has large gravitational 
focusing factors at t = 0, each large planetesimal rapidly sweeps up the small planetesimals along 
its orbit. With growth times proportional to the orbital period (eq. |5-llaj ). protoplanets at 3-7 
AU grow much faster than those at larger a. Growth produces many isolated mass objects packed 
closely together. Early on, gravitational interactions among these objects jostle them around into 
overlapping orbits. After ~ 0.3 Myr, the most massive of these protoplanets begin to scatter lower 
mass protoplanets to smaller and larger a. Scattered protoplanets sweep up and scatter the large 
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pre-existing planetesimals in these orbits, accelerating the growth of all large protoplanets. At ~ 1 
Myr, the largest protoplanets begin to accrete gas from the disk. As they grow, they scatter lower 
mass protoplanets to larger and larger a; eventually, they eject some of these low mass protoplanets 
from the planetary system. 

At the end of the calculation at 100 Myr, six planets remain on stable orbits. Two gas giants, 
with 1 and 3 Jupiter masses, have a = 5 AU and 10 AU. Inside these gas giants, a super-Earth 
with m ~ 10 M0 lies on a fairly circular orbit, e ~ 0.02 at a ~ 1.5 AU. Outside the gas giants, 
two more super-Earths occupy orbits in a 2:1 resonance. After many exchanges, the orbits of these 
two planets are likely stable. Finally, a planet with roughly 1.5 times the mass of Saturn rests in 
an orbit with modest eccentricity, e w 0.1 at a ~ 50 AU. The outcome of this simulation combines 
some properties of the Solar System - four gas giants at 5-30 AU - along with some properties of 
known exoplanets - a super-Earth at 1-2 AU. 

This example illustrates several important differences between terrestrial and gas giant planet 
formation. 

• Oligarchs form at 1 AU before they form at 5 AU. From eq. ( j5-llaj ). the growth time 
scales with the orbital period and the enhancement of the surface density at the snow line. 
Comparing timescales at 0.4 AU and at 4 AU, we expect a factor of ~ 300 from the orbital 
period and a factor of 1/2.4 from the snow line enhancement. The factor of ~ 10 ratio of 
growth times at 0.4 AU and 4 AU in the simulations agrees well with expectations. 

• Oligarchs reach chaotic growth faster at 5 AU than at 1 AU. From the discussion in §5.1.2| 
icy oligarchs at 5 AU have larger isolation masses than rocky oligarchs at 1 AU (eq. |5-12| ). 
Thus, their gravitational interactions are stronger and lead to chaotic growth sooner. 

• While rocky planets scatter low mass protoplanets a few tenths of an AU, gas giant planets 
scatter some low mass protoplanets close to the host star and eject many others from the 
planetary system. Fortunately, the Solar System avoided either outcome. However, many 
exoplanets close to their host stars have large e. Although orbital migration probably accounts 
for exoplanets with small a and nearly circular orbits (Chapter by Morbidelli) , producing 
super-Ea rths with a < 0.4 AU and e > 0.1 probably re quires planet-planet scattering, as in 
Fig. [TT] (j.Turic fc Tremainel lioosl : lEavmond et al.ll201lh . In a few years, exoplanet statistics 



will allow critical tests of the ability of migration and scattering to explain the observed (a, e) 
distribution. 



5.3. Direct Formation of Brown Dw^arfs and Gas Giants 



Though the core accretion is remarkably successful at explaining the diversity of planetary 
systems, it cannot be proved. It is thus useful to develop alternate theories. Ideally each theory 
would make testable predictions, but consistency with physics and existing astronomical observa- 
tions is also provides stringent constraints. Here we consider the leading alternate theory, that 
a gravit ationally unstable gas disk might fragment into bound obiects t hat survive as gas giant 
planets ( Kuiper 1951 : Cameron 1978 : Bodenheimer et al. 1980; Boss 2000l ). This theory should not 
be confused with the hypoth esis that the solid component of the proplanetary disk gravit ationally 
collapses into planetesimals (jSafronovlll969l : lOoldreich fc Ward! [l973l : I Youdin fc Shull20oi Ulh . 
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Section [3T6l introduced the idea that disks are gravitationally unstable when the lToomrd (j 19641 ) 
criterion, 

Q = ^<1, (5-40) 

is satisfied for a surface mass density, H. There, we outlined basic constraints on the accretion 
rate and temperature for stable disks. For a fragment to survive it must cool quickly, so that 
it contracts on an orbital timescale. Anal ytic theory suggests bound fragmen ts are possible but 
have the typical masses of b rown dwarfs (IRafikovlboOSl: Ixratter et al.1 lioiObl) : so far, numerical 
simulations are inconclusive (jD'Angelo et al.ll20ld : ICai et al.ll20ld ). 

To outline the basic issues for gravitational instability, we consider a density perturbation in 
a disk which satisfies the Toomre Q criterion. For this fragment to continue to contract, it must 
cool on a sufficiently short timescale ( Gammi3l200ll ). 



tc < ^ft' 



(5-41) 



For likely cond itions within a protoplanetary disk, numerical simulations suggest that the critical 
value of ~ 3 ( Gammie 2001 : Rice et al. 20031 ) . with uncertainties due to the equation of state and 
the opacity. Fragments that cannot cool on the orbital timescale will be sheared apart. For the 
cooling time in a viscous disk with an a- viscosity (eq. |3-14| ). this constraint places a limit on a: 



a 



> 



4/[9C7(7 - 1)] ~ 0.3. 



(5-42) 



Fragmenting disks also have very high accretion rates. 

For the inner disk (R < 10 AU), combinirig the Toomre (Q < 1) and c ooling criteria yields 
fragments only in hot, massive disks (jRafikovl l2005l : iMatzner &: LevinI l2005l ) . With disk masses 
comparable to the stellar mass, these conditions probably produce a bound stellar companion 
instead of a lower mass planet. Although colder, lower mass disks can have Q ^ 1, their slow 
cooling times do not allow the fragment to survive. 

Forming planet-mass fragments in the outer disk, at 50 - 100 AU, is more attractive. Irradiatio n 
(instead of viscous transport) then dominates the energy budget of the disk ( D'Alessio et al. 1 119981 ). 



Although estimates for the cooling time are more complicated, an irradiated disk is thoug ht to 
fragment more readily once it reaches Q <1 (|Rice et al.ll201ll : iKratter fc Murrav-Clavll201lh . 



In any part of the disk, the initial mass of a gravitational fragment is typically ~ ^gHg in 
terms of the disk scale-height Hg. For the most optimistic assumption s about coolin g - specifically 
that the optical depth r = 1 - fragment masses are > 5 Mj at 100 AU (jRafikovll2005l ). Less efficient 
cooling (r 7^ 1) and closer orbital separations increase the fragment mass. These minimum masses 
require very cold disks (T < 10 K), which might be achievable i n the outer portions of disks in low 
mass star forming regions like Ophiuchus ( Andrews et al. 20091 ). 

While making a fragment with a mass below ~ 10 Mj is challenging, keeping the fragment 
mass low is an even more difficult problem. Disks are likely to fragment before infall from the 
surrounding molecular cloud ceases. Preventing the disk fragment from accreting this infalling 
mass is a challenge. As explained in ^5.2.41 stopping the flow of disk gas onto a planet requires a 
low mass and fairly q uiescent disk, exactly opposite to the conditions required for an unstable disk. 
Kratter et al.l ( 2010bl ) quantifies these issues and concludes that disks in nearby star- forming regions 



are more likely to produce 25-75 Mj brown dwarfs than 1-10 Mj planets. The mass problem might 
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be helped if fragments cool just rapidly enough to remain boun d. Then if they rn i grate inward the y 
would overflow their Roche lobes and be "tidally downsized" ( Bolev et al. 20 id : Navakshin 20101 ) . 
This intriguing suggestion cannot yet be considered a solution. 

Even if bound Jupiter-mass fragments form, understanding the time evolution of fragments 
within a Q ~ 1 disk requires sophisticated numerical calculations. In smooth particle hydrodynam- 
ics (SPH), an ensemble of particles represents the gas; each par ticle has a set of physical propertie s 
and responds to gravity, radiation pressure, and other forces ( Benz et al. 19861 : Monaghan 1992 ). 
Grid-based calculations lay out a set of points where the physical conditions of the gas are specified; 
solvin g a set of coupled differential equations yields the time evolu tion of the conditions at each 
point (|Black k Bodenheimejll975l : lTohlinelE98ol : IPickett et aT]ll998l ). 



In both approaches, d isks close to the stability limit develop multi-arm spiral struc ture (IKr after et al 



2010al: iMaver et al.l 20021 ). Spiral modes generate turbulence throughout the disk (iNelson et al 



19981 : iGammid 200lh. The large a mplitudes of these modes create a rippled surface above the 



disk midplane ( Durisen et al. 200ll ). When the cooling rate is near the critical value ^, the spiral 
structure maintains a rough equilibrium, where the arn plitude of the spiral modes increases as the 
coohng rate dechnes ( Mejia et al. 20051 : Cai et al. 20061 ). 

When the c ooling time is s maller than the local rotat i onal period, as i n equation (j5-4ip . the 



disk fragments ( Gammie 2001 : Johnson Sz Gammie 2003 : Rice et al. 20031 ). For astrophysically 



relevant conditions, fragmentation requires the disk mass be at least 10% of the stellar mass. As 
the disk mass grows, spiral waves propagate through out the disk. Wit hin a few rotation periods, 
fragments form in the densest portion^ of the spirals feolev et a1.I B . If the fragments continue 
to cool rapidly, th ey grow in mass and become bound objects, otherwise they are sheared apart 



(|Maver et al.ll2004l ). 



The long-term evolution of bound fragments in self-gravitating disks is unclear. Because the 
numerical calculations involve such a wide range of scales, none can evolve a bound fragment 
long enough to determine its final fate. If the cooling time is short, and accretion from the disk 
inefficient bound fragment could become gas giant plane ts, or more likely bro wn dwarfs. Such 
objects might also grow to full-fledged stellar companions (|Bonnell fc Batelll994l ). Overcoming the 
numerical limits on the calculations requires faster computers and innovative techniques to follow 
the evolution of fragments in a time-dependent disk. 



6. SUMMARY 



In the last decade or two, observations have revolutionized our understanding of planetary 
astronomy. In the 1990's, the thril ling discovery of th e Kuiper belt nearl y doubled the einpirica l 
size of the disk of the Solar System ( Jewitt Sz Luu 19931 ). A few years later. Mayor Sz Queloz (|l99,5l ) 
discovered an extraordinary exoplanet orbiting only 0.05 AU from the solar-type star 51 Peg. Now, 
the number of known Kuiper belt objects easily exceeds 1000, including a grand variety of dynamical 
and taxonomic cla sses that place inter esting constraints on the origin and early evolution of the 



protosolar nebula ( Barucci et al. 20081 ) . Although there are only ~ 1000 confirmed exoplanets. 



data from Kepler and many ground-based programs will certainly push the count past 10,000 in 
the next decade. Some of these will certainly challenge current ideas about planet formation. 

Throughout this onslaught, theorists responded quickly with new ideas (dead zones in disks, 
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migration, symplectic A^-body codes) and variants on old ideas (collisional cascades, disk instabil- 
ities, multiannulus coagulation codes). Rapid developments in computing hardware fueled many 
advances; new analytical approaches drove others. 

Today, we have a good basic theory of disk evolution. Despite uncertainties about the initial 
mass and temperature distribution and the origin of disk viscosity, analytic and numerical disk 
models provide a good framework for interpreting observations and for exploring the origin of 
planetary systems. Current research involves combining more elaborate versions of the basic theory 

with detailed models for the chemical evolution of disk material. Within the next decade, these 
investigations should improve our insight into the overall evolution of the disk and the growth and 
composition of dust grains with sizes ~ 1-10 mm. 

We also understand how km-sized or larger planetesimals become planets (^. Although there 
are major uncertainties about the onset and the end of gas accretion and the interactions of massive 
planets with the gas disk, analytic theory and numerical simulations demonstrate that - on 1-10 
Myr timescales - ensembles of planetesimals can evolve into terrestrial and gas giant planets within 
~ 50 AU of the central star. Comparisons of observations with the predicted masses and orbital 
properties of planets and the predicted dust masses and luminosity evolution of debris disks are 
promising and will eventually produce stringent tests of the theory. 

Despite these successes, we are still in search of a robust theory for planetesimal formation (SI- 
Excellent progress on the meter-size barrier isolates the importance of radial drift and the problems 
associated with direct coagulation models. Some type of instability - either by direct gravitational 
collapse or a concentration mechanism such as the streaming instability - seems necessary to 
produce planetesimals on fast timescales. Larger numerical simulations will undoubtedly yield a 
better understanding of these instabilities. Exploring other physical mechanisms for particle growth 
and evolution is also essential. 

Understanding fragmentation in protostellar disks promises to unify our understanding of star 
and planet formation. Although many physical and numerical issues remain unresolved, fragmen- 
tation is a promising way to produce gas giants and brown dwarfs at > 50-100 AU from the parent 
star. Despite the current lack of large samples of planets in t his domain, direct imaging surveys are 
starting to discover 1 -30 Mj objects w i th a ^ 10-100 AU (jXalas et allboosl : [ivlarois et ahlbooi : 



Lagrange et~aIll2O10l : ICurrie et al]l2010l : iKraus et allboilh . With large uncertainties in model at- 



mospheres, assessing the formation mechanism of these objects is difficult. Once large samples of 
planets and brown dwarfs at 10-100 AU are available, their properties will allow a robust assessment 
of the core accretion and disk instability mechanisms. 
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Fig. 1. — Schematic view of two adjacent annuli in a disk surrounding a star (black point at left). 
Annulus 1 lies inside annulus 2; material in annulus 1 orbits the star more rapidly than material 
in annulus 2 (r^i > ^2)- 
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Fig. 2. — Time evolution of the surface density for a ring with constant viscosity (eq. |3-5j ). Over 
time, viscous diffusion spreads the ring into a disk. The legend indicates the scaled time, r, for 
each curve. 
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Fig. 3. — Time evolution of the surface density o f a gaseous disk s urrounding a 1 Mq star. Dashed 
hnes show results for the analytic disk model of [Chambers ( 20091 ) : solid lines show results for our 
numerical solution of the diffusion equation. Depsite small differences in the initial conditions, the 
numerical solution tracks the analytic model. 
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Fig. 4. — Time evolution of the disk mass (upper panel) and disk accretion rate onto the central 
star (lower panel) for the analytic and numerical solutions in Fig. |3l 
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Fig. 5. — Aerodynamic stopping time normalized to the Keplerian orbital frequency for a range of 
particle sizes in our reference minimum mass disk model. Small (large) values of Tg indicate strong 
(weak) coupling of solids to the gas disk. The breaks in the curves are due to transitions between 



different drag laws, as described by eq. (j4-3p . An internal density of p, 
the solids. 



1 g cm is assumed for 
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R[/K\J] 



Fig. 6. — Radial drift timescales, R/R, for the same disk models and particle sizes as in Fig. [H 
The fastest drift timescale is for the particle size that has Ts = 1 as indicated by the grey curve. 
Drift timescales are much faster than disk lifetimes of a few Myr, especially near the "meter-sized" 
barrier. 
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Fig. 7. — Disruption energy, Q^, for icy objects. The solid curves plot typical results de- 
rived from numerical sirn ulations of collisions (e.g., Benz Asphaug 19991 : Durda et al. 2004 : 
Leinhardt fc Stewart! [2OO9I ) that include a detailed equation of state for basalt (rock) or crystalline 
ice (ice). In the strength regime (r < lO'^-lO^ cm), smaller particles are stronger. In the gravity 
regime (r > 10^ cm), larger objects are s tronger. The "Rubble P ile" curve shows results consistent 
with model fits to comet breakups (e.g., Asphaug &: Benz 19961 ). 
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Fig. 8. — Evolution of oligarchs in the terrestrial zone. The calculation starts with 1 km plan- 
etesimals (p, = 3 g cm^^) in a disk with Eg = 8 g cm~^ (a/1 AU)~^. Left panel: The time 
evolution of semimajor axis shows three phases that start at the inner edge of the grid and propa- 
gate outward: (i) after runaway growth, isolated oligarchs with m > 4 x lO'^^ g form and continue 
to grow very rapidly; (ii) oligarchs develop eccentric orbits, collide, and merge; and (iii) a few 
massive oligarchs eventually contain most of the mass and develop roughly circular orbits. The 
legend indicates masses (in M^) for the largest oligarchs. Right panel: The mass evolution of 
oligarchs shows an early phase of runaway growth (steep tracks) and a longer phase of oligarchic 
growth (flatter tracks), which culminates in a chaotic phase where oligarchs grow by captures of 
other oligarchs (steps in tracks). Despite the steeper appearance of some of the mass tracks during 
runaway growth, protoplanets grow more rapidly during oligarchic growth. 
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Fig. 9. — Evolution of debris disks in the terrestrial zone (|Kenvon fc BromlevI lionibl . bnOfil ). For 
an A- type star with a luminosity of ~ 50 Lq, the range in blackbody temperatures of planetesimals 
at 3-20 AU (425-165 K) is similar to the range in the Solar System at 0.4-2 AU (440-200 K). 
Left panel: Images of a disk extending from 3-20 AU around an A-type star. The intensity scale 
indicates the surface brightness of dust, with black the lowest intensity and white the highest 
intensity. Right panel: Mid-IR excess for two debris disk models. The light grey line plots the ratio 
of the 24 /im flux from a debris disk at 0.4-2 AU disk relative to the mid-IR flux from a G-type 
star. The dark grey line shows the evolution for the A-star disk shown in the left panel. 
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Fig. 10. — Mass of planet atmosphere as a function of core mass and accretion time (eq. |5-36| ). 
The legend indicates the accretion time in Myr. At fixed core mass, the mass of the atmosphere 
grows with the accretion time. Longer accretion times allow more massive cores to have hydrostatic 
atmospheres. 
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Fig. 11. — Orbital evolution of icy oligarchs. The calculation starts with 1 cm and 1000 km 
planetesimals {p, = 1.5 g cm~^) in a disk with initial Eg = 14g cm~^ AU)~"'^e~"/'^'^ . 

During the first 3 x 10^ yr, ~ 20 oligarchs with m ~ 0.05 form in a relatively narrow range of 
semimajor axis, 3-7 AU. As they grow, more massive oligarchs scatter lower mass oligarchs to large 
semimajor axes, a ~ 1-20 AU. At ~ 1 Myr, some oligarchs begin to accrete gas. Over the next 
10 Myr, continued growth and scattering leads to collisions and mergers of oligarchs. Eventually, 
only a few oligarchs remain. The largest of these have masses comparable to the mass of Jupiter. 
The legend indicates masses (in M^) for the largest oligarchs in stable orbits around the central 
star. 
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Table 1:: Frequently Used Symbols 



Symbol Ref. Meaning 



General Physical Quantities 

AU eq. ()2-ip astronomical unit 

c eq. ()2-8p speed of light 

Cs eq. ()3-3p sound speed 

Cd eq. ()4-4p drag coefficient 

Cy H'S.2\ specific heat 

P eq. (|3-2p gas pressure 

t time or timescale, often subscripted 

T [Tcq] eq. (|2-5p [equilibrium] temperature 

vk eq. (jS-Sp Keplerian circular velocity 

7 eq. (|3-3p adiabatic index 

K eq. (|3-33p radiative opacity 

A eq. ()2-6p . ^3.61 wavelength (light or other waves) 

A §3.21 eq. (I4-3P gas mean free path 

fj, H3.2\ mean molecular weight 

a ^3.21 cross section 

usB eq. ()3-23|) Stephan-Boltzmann constant 

p density (mass per volume) of quantity, often subscripted 

p, eq. (l2-8p Internal density of solid grain or planetesimal 

f2 ^3.11 Keplerian orbital frequency 



Stellar Quantities 

Lj, [Lq] eq. (l2-5p stellar [solar] luminosity 

l-^: eq. (|5-26p L^/Lq, dimensionless stellar luminosity 

[Mq] eq. (f2^ stellar [solar] mass 

Mj,/M0, dimensionless stellar mass 

eq. (12-9P stellar metallicity 

Disk Quantities 

D eq. (j3-18p Dissipation (i.e. heating) rate, per area 

F eq. (f2^ Disk mass relative to MMSN 

if eq. (|3-13p disk scaleheight 

J [J] eq. (j3-10p angular momentum [torque] 

Luce [Ld] eq. (l3-2ip [eq. (13-19p ] Accretion luminosity [part released in the disk] 

M eq. ()3-ip accretion rate through disk 

Q eq. ()3-34p dimensionless measure of gravitational stability 

R eq. (|2-ip distance from star (cylindrical radius) 

i?in eq. (I3-19P radius of disk inner edge 

R J eq. (|3-17p torque on inner disk edge, as length-scale 

v,f, [vr] eq. (13- ip orbital [radial] flow velocity 

Zfjisk eq. (|2-2p Disk "metallicity" as ratio of solids to gas 
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Table 1:: Frequently Used Symbols 



Symbol 


Ref. 


Meaning 




eq. (12^41) 


"Metallicity" relative to fiducial Solar value of 0.015 


a 


eq. (|3-14|) 


dimensionless angular momentum transport coefficient 


OtD 


eq. (|4-12|) 


H 1 TTi f^Ti m nn 1 i"!! y ni 1 1 pttI' n itti n t;i on rri ptti r*! pnf 


V 


eq. (FTD 


fraction by which rotation is slower than Keplerian 


e 


eq. (I3-26D 


grazing angle of starlight on disk surface 


V 


eq. (p^^I) 


viscosity, usually "anomalous" 


PR 


eq. (j4-13|) 


Roche density for gravitational binding 




eq. (12^. eq. (IH^T]) 


Surface density (mass per area) of disk, e.g. gas (§3]) or plan- 






etesimals (^5.ip. Subscripted as needed. 




eq. (M [eq. (l2^] 


Surface density of gas [particle! disks. 


t 


eq. (la^ij) 


Inflow or outflow of mass from disk, per area 


Protoplanet Quantities 




a 


eq. (IFTl) 


semimajor axis (similar to disk R) 


B 


eq. (|5-12j) 


width of feeding zone in Rh 


fc 


ea. (IF7D & (15-lUp 


gravitational focusing factor for collisional cross section 


m [nic , ma] 


gSjJ [^5.2.11 


protoplanet mass, total or [core, atmosphere] 


rus [mi] 


m\ 


mass of small protoplanets, i.e. planetesimals [larger proto- 






planets] 


miso 


eq. (I5-12D 


isolation mass 


m [rh;] 


§5.1.21 


accretion rate, for [large] protoplanet's mass growth 


Q*D 


eq. (|5-18D 


energy (per mass) threshold for catastrophic collisional dis- 






ruption 


r [ts , ri , rc] 


gSD, §5.2.11 


radius of protoplanet [small, large, solid core] 




eq. (|5-24[) 


Bondi radius for planet gravity to exceed thermal energy of 






disk gas 


Rh 


eq. (|FT1) 


Hill radius for planet gravity to dominate stellar tides 


s 


eq. (li^Hj) 


radius of small grain or planetesimal 


ts [Ts] 


ea. (13:31) [ea. (H^SD] 


[dimensionless] aerodynamic stopping timescale 


V [Vs , Vl] 


ED 


velocity dispersion of [small, large] protoplanets 


Vesc [vh] 


ea. ([Fi]) [ea. (IF3])] 


escape speed from protoplanet's surface 


[I^l] 




surface mass density of small [large] protoplanets 


X 


eq. (|5-:-i5[) 


log ratio of atmosphere to core radius 






r/Rff, radius in Hill units 
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